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The  following  conditions  regulate  the  use  of  computer  programs 
developed  by  the  Hydrologic  Engineering  Center  (HEC),  Corps  of 
Engineers,  Department  of  the  Army. 

1.  The  computer  programs  are  furnished  by  the  Government  and  are 
accepted  and  used  by  the  recipient  individual  or  group  entity 
with  the  express  understanding  that  the  United  States  Government 
makes  no  warranties,  expressed  or  implied,  concerning  the 
accuracy,  completeness,  reliability,  usability,  or  suitability 
for  any  particular  purpose  of  the  information  or  data  contained 
in  the  programs,  or  furnished  in  connection  therewith,  and  that 
the  United  States  Government  shall  be  under  no  liability 
whatsoever  to  any  individual  or  group  entity  by  reason  of  any  use 
made  thereof . 

2.  The  programs  belong  to  the  United  States  Government. 
Therefore,  the  recipient  agrees  not  to  assert  any  proprietary 
rights  thereto  nor  to  represent  the  programs  to  anyone  as  other 
than  Government  programs . 

3.  The  recipient  may  impose  fees  on  clients  only  for  ordinary 
charges  for  applying  and  modifying  these  programs. 

4.  Should  the  recipient  make  any  modifications  to  tne 
program(s),  the  HEC  must  be  informed  as  to  the  nature  and  extent 
of  those  modifications.  Recipients  who  modify  HEC  computer 
programs  assume  all  responsibility  for  problems  arising  from,  or 
related  to,  those  modifications.  User  support  from  the  HEC  to 
third  party  recipients  will  only  be  provided  after  the  second 
party  demonstrates  that  program  difficulties  were  not  caused  by 
their  modifications. 

5.  This  "Conditions  of  Use"  statement  shall  be  furnished  to  all 
third  parties  that  receive  copies  of  HEC  programs  from  the 
recipient.  Third  party  recipients  must  be  notified  that  they 
will  not  receive  routine  program  updates,  correction  notices,  and 
other  program  services  from  the  HEC  unless  they  obtain  the 
program(s)  directly  from  the  HEC. 

6.  All  documents  and  reports  conveying  information  obtained  as  a 
result  of  the  use  of  the  program(s)  by  the  recipient,  or  others, 
will  acknowledge  the  Hydrologic  Engineering  Center,  Corps  of 
Engineers,  Department  of  the  Army,  as  the  origin  of  the 
program ( s ) . 
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1.  INTRODUCTION 

This  program  was  prepared  in  the  Hydrologic  Engineering  Center. 
Up-to-date  information  and  copies  of  source  statement  cards  for 
various  types  of  computers  can  be  obtained  from  the  Center  upon 
request  by  Government  and  cooperating  agencies.  While  every  care 
is  taken  to  validate  this  program,  it  is  not  feasible  to  anticipate 
and  test  all  possible  applications.  Consequently,  the  Center  is 
interested  in  problems  that  arise  in  application  and  will  assist 
in  resolving  deficiencies  in  the  program  to  the  extent  feasible. 

2.  PURPOSE  OF  PROGRAM 

The  purpose  of  this  program  is  to  perform  frequency  computations 
of  annual  maximum  hydrologic  events  necessary  to  a  regional  frequency 
study.  Frequency  statistics  are  computed  for  recorded  events  at  each 
station  and  for  each  duration.  Missing  events  are  computed  so  that 
complete  sets  of  events  are  obtained  for  all  years  at  all  stations 
while  preserving  all  inter-correlations.  These  are  arranged  in  the 
order  of  magnitude  for  each  station  and  duration  and  tabulated  with 
median  plotting  positions.  Statistics  for  each  station  are  then 
adjusted  to  the  complete  period  of  region  record,  and  frequency 
curves  are  computed  in  accordance  with  procedures  given  in  "Statistical 
Methods  in  Hydrology"  by  Leo  R.  Beard,  January  1962 ,  using  the  logarith¬ 
mic  Pearson  Type  III  function  and  the  expected -probability  concept. 

The  use  of  all  long-record  stations  instead  of  only  one  for  the  extension 
of  frequency  statistics  at  short-record  stations  is  considered  to  consti¬ 
tute  some  advantage  over  procedures  given  in  "Statistical  Methods" .  As 
an  alternative  use  of  this  program,  frequency  statistics  can  be  supplied 
and  curves  will  be  computed. 

3.  DESCRIPTION  OF  EQUIIMENT 

A  FORTRAN  IV  compiler,  random  number  generator  (function  RNGEN 
included,  see  Exhibit  2),  and  large  memory  are  required.  The  large 
amounts  of  computation  make  high  speed  desirable.  Accordingly,  it  is 
virtually  necessary  to  use  a  computer  of  the  IBM  709U  class  for 
execution  of  this  program.  It  is  desirable  to  use  one  input  tape  and 
one  output  tape  unit,  in  addition  to  card  (tape  7)  and  printer  (tape 
6)  output  and  standard  (tape  5)  input. 


4.  METHODS  OF  COMPUTATION 


a.  Flows  for  those  stations  with  zeros  in  the  data  are  first  incremented 
by  l/lO  percent  of  their  average  for  each  station  and  duration  in  order  to 
preclude  infinite  negative  logarithms.  This  increment,  if  added,  is  later 
subtracted  from  reconstituted  flows  and  computed  frequency  curves.  The  mean, 
standard  deviation  and  skew  coefficient  of  the  logarithms  for  each  station 
and  duration  are  then  computed.  Preliminary  to  estimating  missing  flows 
by  correlation,  each  flow  is  then  converted  to  a  standardized  variate  using 
an  approximation  of  the  Pearson  Type  III  distribution.  This  involves  the 
following  equations: 
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in  which: 

X  =  Logarithm  of  flow  event 
Q  =>  Recorded  flow  event 

q  =  Small  increment  of  flow  used  to  prevent  infinite 
_  logarithms  for  events  with  zero  flow 
X  =  Mean  logarithm  of  flow  events 
N  a  Total  years  of  record 

S  =  Unbiased  estimate  of  population  standard  deviation 
g  =  Unbiased  estimate  of  population  skew  coefficient 
t  =  Pearson  Type  III  standard  deviate 
i  a  Duration  number 
m  =  Year  number 
K  =  Normal  standard  deviate 
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b.  After  transforming  the  flows  for  all  stations  and  durations 
to  normal,  the  gross  (simple)  correlation  coefficients  R  between  all 
pairs  of  stations  for  each  duration  and  for  adjacent  durations  at  each 
station  are  computed  by  use  of  the  following  formula: 

f  r  N  o  N  o  w  o 

R  =  -  [l  -  (  2  (x  x  )2/(2  x2  _  2  X?  ) 

i  1  t  m-l  i,m  x-l,m  m=l  i,m  m-1  x-l,m  J 


(N-l)/(N-2) 


(7) 


c.  Inasmuch  as  not  all  stations  and  durations  necessarily  have 
the  same  length  of  record,  correlation  matrices  obtained  in  b  might 
not  be  complete  or  internally  consistent.  If  not,  missing  values 
are  estimated,  and  low  values  are  raised  to  obtain  consistency, 
inasmuch  as  lew  values  are  least  reliable  and  least  influential. 

Each  missing  value  is  estimated  by  examining  its  relationship  to 
related  pairs  of  values  by  use  of  the  following  formula,  using  i , 
j  and  k  subscripts  to  indicate  variables  used  in  the  gross  correlation: 


Rij  RkiRkj 


2)(l-RkJ2) 


(8) 


d.  Consistency  of  each  correlation  matrix  to  be  used  for 
estimating  missing  flows  is  assured  by  first  testing  all  combinations 
of  triads  of  correlation  coefficients  used  in  that  matrix.  The  test 
for  consistency  of  each  complete  matrix  is  made  by  computing  the 
multiple  con-elation  coefficient.  If  this  value  is  greater  than  1.0, 
further  adjustment  is  required.  Such  further  adjustment  is  obtained 
by  introducing  a  coefficient,  successively  smaller  by  .2,  on  the 
radical  in  equation  8  and  repeating  all  triad  consistency  tests  until 
all  matrices  are  consistent. 

e.  Missing  flows  are  estimated  by  correlation  with  corresponding 
flows  at  other  stations  and  the  flow  at  the  same  station  for  the  adjacent 
duration  (preceding  duration,  except  that  the  succeeding  duration  is 
used  when  estimating  for  the  first  tabulated  duration).  Since  it  is 

not  known  which  stations  might  have  recorded  or  previously  estimated 
values,  the  correlation  matrix  and  regression  equation  might  be  different 
for  the  same  station  and  duration  in  different  years.  The  regression 
equation  is  computed  for  each  missing  value  in  terms  of  normal  standard 
variates  by  selecting  required  coefficients  from  the  complete  (and 
consistent)  correlation  matrix  and  solving  by  the  Crout  method  explained 
in  Exhibit  1.  The  missing  value  Is  computed  from  this  regression  equation. 
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introducing  a  random  component  equal  to  the  non-determination  of 
the  equation,  in  order  to  preserve  the  proper  variance  (standard 
deviation)  of  the  flows.  This  is  done  as  follows: 


*1  =  @2^2  +  ^3k3  + 


(9) 


in  which: 

k  =  Normal  standard  deviate 

/3  =  Beta  coefficient 

FT  =  Determination  coefficient 

Z  =  Random  number  normally  distributed 

n  =  Number  of  variables  in  equation 

f.  When  all  flows  have  been  reconstituted,  the  mean  and  standard 
deviation  for  each  station  and  duration  are  recomputed.  Regression 
lines  of  standard  deviation  and  skew  coefficient  separately  versus  mean 
are  computed,  and  "smoothed"  values  of  standard  deviation  and  skew 
obtained  as  described  in  "Statistical  Methods" .  Equivalent  record  for 
the  recorded  and  reconstituted  flows  for  each  station  and  duration  is 
estimated  by  adding  the  determination  coefficient  for  each  year  of 
reconstituted  flow  to  the  total  years  of  recorded  flews.  This  equivalent 
record  is  used  in  computing  expected  probabilities  as  discussed  below. 
Flows  are  arranged  in  descending  order  of  magnitude  and  median  plotting 
positions  are  computed  as  defined  in  "Statistical  Methods".  Frequency- 
curve  coordinates  for  each  station  and  duration  are  computed  from  the 
mean,  standard  deviation,  skew  coefficient,  flow  increment  and  equivalent 
record  length,  using  table  values  of  the  normal  distribution,  the  trans¬ 
form  for  the  Pearson  Type  III  function  shown  in  Equation  10,  and  the 
following  approximate  transforms  for  expected  probability: 

P  „  =  .01  (1+1600/N1*72)  (10) 

•Ux 

P  =  .1  ( 1+280 /N1,55)  (11) 

•JL 

p^^  =  1  +  26/N1,16  (12) 


u 


in  which: 


'5  =  5  (1  -!  6/N1*04) 

(13) 

>1Q  =  10  (l  +  3/N1*04) 

(lk) 

>3Q  =  30  (1  +  .46/N*925) 

(15) 

P  =  Expected  probability  in  percent,  symmetrical  about  50  percent 
N  *s  Equivalent  years  of  record 

5.  INPUT 

Input  is  suiranarized  in  Exhibits  6  and  J.  All  data  are  entered 
consecutively  on  each  card,  using  8  columns  (digits,  including 
decimal  point,  if  used)  per  variable  and  10  variables  per  card  unless 
fewer  variables  are  called  for,  except  that  the  first  column  on  each 
card  is  reserved  for  identification,  ftie  first  output  title  card 
must  have  an  A  in  column  1.  An  example  of  input  is  given  in  Exhibit 
2.  Certain  inadequacies  of  data  will  abort  the  job  and  waste  input 
cards  until  the  next  card  with  A  in  column  1  is  reached.  After  a 
job  is  finished,  a  card  with  A  in  column  1  followed  by  3  blank  cards 
causes  the  computer  to  stop. 

6.  OUTPUT 

Printed  output  includes  key  input  information  for  job  identification 
and  all  results  of  computations.  An  example  of  printed  output  is  given 
in  Exhibit  3. 

7.  OPERATING  INSTRUCTIONS 

Standard  FORTRAN  IV  instructions  and  random  number  generator  are 
required.  No  sense  switches  are  used. 

8.  DEFINITIONS  OF  TERMS 

Terms  used  in  the  program  are  defined  in  Exhibit  4. 

9.  PROPOSED  FUTURE  DEVELOPMENT 

No  specific  future  development  of  this  program  is  presently  planned. 
It  is  requested  that  any  user  who  finds  an  inadequacy  or  desirable  addi¬ 
tion  or  modification  notify  the  Hydrologic  Engineering  Center. 
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EXHIBIT  1 
Crout's  Method 


One  of  the  best  methods  for  solving  systems  of  linear  equations 
on  desk  calculating  machines  was  developed  bv  P.  D.  Crout  in  1941. 

This  method  is  based  on  the  elimination  method,  with  the  calculations 
arranged  in  systematic  order  so  as  to  facilitate  their  accomplishment 
on  a  desk  calculator.  In  this  method  the  coefficients  and  constant 
terms  of  the  equations  are  written  in  the  form  of  a  "matrix,"  which  is 
a  rectangular  array  of  quantities  arranged  in  rows  and  columns. 

The  method  is  best  explained  by  an  example.  Suppose  that  in  a 
multiple  correlation  analysis  it  is  required  to  solve  the  following 
system  of  linear  equations  to  obtain  the  unknown  values  of  b  ,  b,,  b, 
and  b^  5 
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For  simplicity  let  us  replace  the  coefficients  of  the  b’s  by  the  letters 
p,  q,  r  and  s,  and  the  constant  terras  by  the  letter  t,  using  subscripts 
1,  2,  3  and  4  to  denote  the  respective  equations! 
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A  continuous  oheok 

on 

the 

computations  as 

they  progress  may  be  obtained 

by  adding  to  the  matrix  of  the  above  system  a  column  of  u*s,  such  that 
u-p  +  q  +  r+  s  +  t.  The  matrix  and  check  column  are  written  as 
follows < 
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The  elements  p^,  q^»  and  form  the  "principal  diagonal"  of 

the  matrix.  Examination  of  the  original  equations  shows  that  the  co¬ 
efficients  are  symmetrical  about  the  principal  diagonal,  i.e.,  q^  «  p2, 

rl  "  p3’  r2  *  q3’  81  "  p4*  s2  *  V  8114  s3  “  V 

This  is  characteristic  of  the  system  of  equations  to  be  solved  in  any 
multiple  correlation  analysis.  Because  of  this  symmetry,  the  computa¬ 
tions  are  considerably  simplified.  While  the  Crout  method  may  be  used 
to  solve  any  system  of  linear  equations,  the  computational  steps  given 
here  are  applicable  only  to  those  with  symmetrical  coefficients. 

The  solution  consists  of  two  parts,  viz.,  the  computation  of  a 
"derived  matrix"  and  the  "baok  solution."  Let  the  derived  matrix  be 
denoted  as  follows! 
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EXHIBIT  1 


The  elements  of  the  derived  matrix  are  oompruted  as  follows t 
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The  general  pattern  of  the  above  computations,  which  may  be  applied 
to  a  system  containing  any  number  of  equations,  is  as  followsi 

(1)  The  first  column  of  the  derived  matrix  is  copied  from  the  first 
oolumn  of  the  given  matrix. 

(2)  The  remaining  elements  in  the  first  row  of  the  derived  matrix 
are  oomputed  by  dividing  the  corresponding  elements  in  the  first  row  of 
the  given  matrix  by  the  first  element  in  that  row. 

(3)  After  completing  the  n—  row,  the  remaining  elements  in  the  (n+l)— 
oolumn  are  oomputed.  Suoh  an  element  (X)  equals  the  corresponding  element 
of  the  given  matrix  minus  the  produot  of  the  element  immediately  to  the  left 
of  (X)  by  the  element  immediately  above  the  principal  diagonal  in  tbs  same 
column  as  (X),  minus  the  produot  of  the  second  element  to  the  left  of  (X) 

by  the  second  element  above  the  prinoipal  diagonal  in  the  same  column  &a  (X), 
eto.  After  each  element  below  the  prinoipal  diagonal  is  recorded,  and  while 
that  element  is  still  in  the  caloulator,  it  is  divided  by  the  element  of 
the  prinoipal  diagonal  which  is  in  the  same  column.  The  quotient  is  the 
element  whose  looation  is  symmetrical  to  (X)  with  respect  to  the  prinoipal 
diagonal. 
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EXHIBIT  1 


(4)  When  the  elements  in  the  (n+1)^  column  and  their  symmetrical 
counterparts  have  been  recorded,  the  (n+l)lfe  row  will  be  aomplete  except 
for  the  last  two  elements,  which  are  next  computed.  Such  an  element  (X) 
equals  the  corresponding  element  of  the  given  matrix  minus  the  produot 
of  Idle  element  immediately  above  (X)  by  the  element  immediately  to  the 
left  of  the  principal  diagonal  in  the  sam  row  as  (X),  minus  the  produot 
of  the  second  element  above  (X)  by  the  second  element  to  the  left  of  the 
principal  diagonal  in  the  same  row  as  (X),  eto.,  all  divided  by  the 
element  of  the  principal  diagonal  in  the  same  row  as  (X). 


The  cheok  column  (u)  of  the  derived  matrix  serves  as  a  continuous 
check  on  the  computations  in  that  each  element  in  the  column  equals 
one  plus  the  sum  of  the  elements  in  the  same  row  to  the  right  of  the 
principal  diagonal.  That  is, 

^  *  1  +  +  S^  +  T^ 


1  +  ft 


S2  +  T2 


U, 


1  +  S-  +  T- 
3  3 


U. 


1  +  T, 


This  check  should  be  made  after  completing  each  row* 

The  elements  of  the  derived  matrix  to  the  right  of  the  principal 
diagonal  form  a  system  of  equations  which  may  now  be  used  to  compute 
the  unknown  values  of  b^,  b^,  b^  and  b^  by  successive  substitution* 

This  is  known  as  the  "back  solution."  The  computations  are  as  follows: 
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S2b5  -  RA 
Slb5  “  Rlb4  “  Qlb3 


It  is  very  important  that  the  computations  be  carried  to  a  suffi¬ 
cient  number  of  digits,  both  in  oomputing  the  coefficients  and  oonstant 
terms  of  the  original  equations,  and  in  computing  the  elements  of  the 
derived  matrix.  It  is  possible  for  relatively  small  errors  in  the 
coefficients  and  oonstant  terms  of  the  original  equations  to  result  in 
relatively  large  errors  in  the  oomputed  solutions  of  the  unknowns.  The 

4 

EXHIBIT  1 


greatest  source  of  error  In  computing  the  elements  of  the  derived  matrix 
arises  from  the  loss  of  leading  significant  digits  by  subtraction.  This 
must  be  guarded  against  and  can  be  done  bv  carrying  the  computations  to 
more  figures  than  the  data.  As  a  general  rule,  it  is  recommended  -chat 
the  coefficients  and  constant  terms  of  the  original  equations  be  carried 
to  a  sufficient  number  of  deoimals  to  produce  at  least  five  significant 
digits  in  the  smallest  quantity,  and  that  the  elements  of  the  derived 
matrix  be  carried  to  one  more  decimal  than  this,  but  to  not  less  than 
six  significant  digits. 
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EXHIBIT  2 


RA'vJMM  NUMBER  FUNCTION  ’INGE'.’ 


This  random  number  function  is  for  a  binary  machine  and  the  constants  must  lie 
computed  according  to  the  number  of  hits  in  an  inteeer  "ord.  The  numbers 
generated  are  uniform.lv  distributed  in  the  interval  0  to  1. 

The  function  is  called  from  the  main  program  by  a  statement  similar  to  the 
following : 


A  =  RMGEN  (IX) 

Where  A  is  some  floating  point  variable  name  and  IX  is  some  integer  variable 
name.  fho  argument  name  IX  need  not  be  the  same  in  the  main  program  aid  the 
function.  The  argument  must  he  initialised  to  zero  in  the  main  progra"’ .  The 
location  of  the  initializing  statement  is  important  and  depends  on  the  results 
desired.  If  it  is  desired  to  have  different  sets  of  random  numbers  for  each 
of  several  different  sets  of  computations  (jobs)  that  are  run  sequent iallv  on 
the  same  program,  then  the.  argument  must  he  initialized  at  the  verv  beginning 
of  the  program  and  never  reinitialized.  Tf  it  is  permissible  to  use  the  same 
sequence  of  random  numbers  for  each  job,  the  argument  must  he  initialized  at 
the  beginning  of  each  job.  The  advantage  of  this  latter  option  occurs  when  one 
of  the  jobs  must  be  re--run  for  some  minor  reason  as  the  sane  random  numbers 
T.’fll  be  used  and  the  results  will  be  comparable. 

Three  constants  must  be  computed  bv  the  following  equations: 

Constant  one  (Cl)  =  +  3 

Constant  two  (C2)  =  2  -1 

>5 

Constant  three  (Cl)  =  1 

'There:  9  =  number  of  bits  in  an  integer  word 


The  constants  for  some  of  the  common  computers  are  listed  in  the  following  table 


CTU’UTi'.R  _ 

Cl!  200  Series 
Ci.  400  Series 
Ids'!  360  Series 
I IV 1  7  )40  and 
7000  Series 
CXI VAC  11 OH 
(.DC  bOOO  Series 


si 7b  or 

I  NT  ECU'! 

WORT) 

Cl 

19 

102  7 

.2  3 

4oqo 

31 

6.55  3'' 

35 

2  60 1  4  7 

4 h  16777210 


COT^TAXTS 
_ 2_  _ 

57.4287 

8 388607 

2147433647 

34359738367 


231474n767r)65.3 


C3 

0 . 1907348630-05 
0. 11920920011-06 
0 . 465661287K-09 
0 . 291038  30460-10 


0.  355271367311-14 
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EXAMPLE  INPUT 


regional  FREQUENCY  COMPUTATION 
TEST  DATA 
JULY  1972 


1 

1945 

1 

PEAK 

0. 

32 

1945 

77100 

32 

1946 

206000 

32 

1948 

185000 

32 

1949 

137000 

32 

1950 

99000 

TEST 

DATA 

723-X6- 

L73SC 

MULTIPLE 

ST  AT  ION 

AND  DURATION 

5 

1945 

PEAK 

1-DAY 

3-DAY 

10-DAY 

30-DAY 

32 

1945 

77100 

71200 

62000 

51000 

30830 

32 

1946 

206000 

185000 

134000 

834C0 

51000 

32 

1947 

133000 

133000 

115000 

65300 

43670 

32 

1948 

185000 

167000 

132000 

85600 

44130 

32 

1949 

137000 

122000 

7040C 

66800 

38130 

32 

1950 

99000 

95900 

90000 

64200 

m6100 

35 

1946 

48400 

32500 

24300 

12870 

7493 

3S 

1947 

46000 

32600 

29270 

16020 

9570 

35 

1948 

53400 

40300 

24870 

12980 

6890 

35 

1949 

18600 

14600 

10570 

8090 

5690 

35 

1950 

23600 

20100 

15800 

9840 

6920 

TEST  i 

DATA 

723-X6-I 

L7350 

SAVE  STATIONS  FROM  PREVIOUS 

JOB 

5 

1945 

1 

2 

PEAK 

1-DAY 

3-DAY 

10-OAY 

30-DAY, 

-.2 

-.4 

-.5 

-.6 

-.8 

32 

35 

33 

1945 

5530 

5040 

4100 

3320 

2270 

33 

1946 

13300 

9560 

7700 

4840 

3150 

33 

1947 

10300 

9360 

8530 

4850 

3540 

33 

1948 

1C300 

3840 

6930 

4230 

2790 

33 

1949 

6470 

5400 

4300 

3120 

2330 

TEST  DATA 

723-X6-L7350 

STATISTICS  FURNISHED 

5 

1945 

1 

PEAK 

1-DAY 

3-OAY 

10-DAY 

30-DAY 

-.2 

-.4 

-.5 

-.6 

-.8 

32 

PEAK 

5.123 

.  159 

-.334 

0. 

6.0 

32 

1-OAY 

5.089 

.153 

-.366 

0. 

6.0 

32 

3-DAY 

4.984 

.133 

-.462 

0. 

6.0 

32 

1 0-OAY 

4.835 

.106 

-.599 

0. 

6.0 

32 

30-QmY 

4.621 

.066 

-.795 

0. 

6.0 

35 

PEAK 

4.519 

.196 

-.278 

0. 

5.6 

35 

l-DAY 

4,408 

.177 

-.168 

0. 

6.0 

35 

3-OAY 

4.267 

.153 

-.027 

0. 

6.0 

35 

10-DAY 

4.052 

.117 

.188 

0. 

5.8 

35 

30-DAY 

3.843 

.082 

.398 

0. 

5.9 

A 

A 

A 
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EXAMPLE  output 


«******«#•>•♦***»*****««***♦**** 

JULY  1QT2  723-X6-L2350 

REGICNAu  FREQUENCY  CONFUTATION 
VERSION  DATE  -  AUGUST  21,  la7S 

it***************************** 


regional  frequency  computation 

TEST  DATA 
JULY  1972 


ndur 

IYRA 

I3KEK 

KEEP 

ICOHV 

IPCHQ 

IPCHS 

NSTAT 

N3MTW 

INCAD 

1 

1905 

1 

-0 

-0 

•0 

-0 

-0 

-0 

-0 

regional  skew  coefficients- 

peak 

o. 


FREQUENCY  statistics  OP  recorded  OATA 


STA 

ITEM 

PEAK 

32 

MEAN 

5.120 

STD  DEV 

,i<70 

SKEW 

*.296 

INCRPT 

0. 

YEARS 

5. 

RECORDED  AND  ^CONSTITUTED  DATA 


3T  A 

YEAR 

PEAK 

32 

19G5 

771C0. 

32 

1936 

206000. 

32 

19CA 

lescoo. 

32 

JR.’.R 

137000. 

32 

19'j0 

99000. 

FREQUENCY 

ARRAYS 

station 

32 

NO  PLOT 

PEAK 

t  :?,*<} 

206000, 

c  3'..<-7 

1 6500  0  . 

3  50. CC 

137000. 

a  60.33 

9  9  0  0  C  • 

5  B  7 .06 

77100. 

AODPTEO  FREQUENCY  STATISTICS 

STA 

ITEM 

peak 

32 

mean 

5,120 

STD  OEV 

.180 

SKEn 

0. 

INCRMT 

0. 

COMPUTEO  frequency  curves 

STATION  3? 


PLOT 

Exp  pons 

PEAK 

.01 

1.01 

617037. 

.10 

2.01 

073707, 

1.00 

5.02 

3«5?33. 

5.00 

10.63 

259953, 

10.00 

15.63 

223966, 

3C.00 

33.11 

1655:8. 

50.00 

50.00 

131853. 

70.00 

66.89 

106320. 

90.00 

80.37 

77620. 

95.00 

89,37 

66378. 

99,00 

9C.98 

50263, 

99, qo 

97.59 

36697. 

99.99 

98,99 

28137. 
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EXHIRIT  4 


JULY  1972  723-X6-L235C 
REGIONAL  FREQUENCY  COMPUTATION 
VERSION  DATE  -  AUGUST  21,  1979 


TEST  DATA 
723-X6-L2350 

MULTIPLE  STATION  A|,0  DURATION 

NQUR  I YRA  ISKSW  KEEP  ICONV  IFCHQ  IPCM5  NSTAT  NSMTH 


5  1945 

-0 

-0 

-0 

-0 

-0 

FREQUENCY  STATISTICS 

OF  RECORDED  data 

STA  ITEM 

peak 

t-OAY 

3-OAY 

10-DAY 

30-DAY 

32  Mc»N 

5.123 

S  ,069 

4.904 

4.835 

4.621 

3TD  OEV 

.161 

.154 

.142 

.C83 

.076 

SKEW 

-.388 

-.527 

-.375 

-.266 

-1.088 

I NCRMT 

0. 

0. 

0. 

0. 

0. 

years 

6. 

6. 

6. 

6. 

6. 

35  MEAN 

4.544 

4.420 

4  •  2q4 

4.066 

3.853 

STO  OEV 

.208 

.181 

.181 

.116 

.082 

SKEW 

-.689 

-.721 

-.964 

-.334 

.593 

INCPMT 

0. 

0. 

0. 

c. 

0. 

YEARS 

5. 

5. 

5. 

5. 

5. 

FREQUENCY  STATISTICS 

AFTER 

ADJUSTMENT 

WITH  A  LONS 

TERM  STATION 

STA  ITEM 

PEAK 

i-day 

3-DAY 

10-DAY 

30-DAY 

32  MEAN 

5.123 

5.089 

4.984 

4.835 

0.621 

STO  OEV 

.161 

.154 

.142 

.083 

.076 

SKEW 

-.334 

-  .  366 

-.462 

-.599 

-.795 

I  NCR  MT 

0. 

0. 

0. 

0. 

0. 

EOUIV  YRS 

6.0 

6.0 

6,0 

6.0 

6.0 

35  MEAN 

4.498 

4.379 

4,242 

4,033 

3.838 

STD  OEV 

.227 

.202 

.206 

.133 

.091 

SKEW 

-.734 

-.612 

-.478 

-.269 

-.073 

INCRMT 

0. 

0. 

0. 

0. 

0. 

EQUT V  YRS 

5.3 

5.3 

5.8 

5.0 

5.0 

INCAD 

-0 
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EXHIBIT  4 


CORRELATION  COEFFICIENTS  OF  RECORDED  OATA  FOR 


PEAK  OURATION 


STA 

32 

35 

WITH  SANE  DURATION 

32 

1.000 

.616 

35 

.M6 

1.000 

WITH  ADJACENT  DURATION 

AT  ABOVE 

STATION 

32 

.995 

.090 

35 

.Tia 

.082 

CORRELATION 

COEFFI 

CIENTS  OF  RECORDED  DATA 

for 

1-OAV  duration 

STA 

32 

35 

WITH  SANE  DURATION 

32 

l.OCO 

.600 

35 

,bO* 

1  .000 

WITH  ADJACENT  DURATION 

AT 

above 

STATION 

32 

.995 

.710 

35 

.420 

.052 

CORRELATION 

COEFFICIENTS  OF  RECORDED  DATA 

FOR 

3-DAV  DURATION 

STS 

32 

35 

WITH  SAMS  DURATION 

32 

1.000 

.867 

35 

.as; 

1.000 

WITH  ADJACENT  DURATION 

AT 

ABOVE 

STATION 

32 

.403 

.909 

35 

,330 

.  893 

CORRS 

L  AT  ION 

COEFFICIENTS  OF  RECOROEO  DATA 

FOR  10- 

DAY  DURATION 

STS 

32 

35 

WITH  SANE  DURATION 

32 

l.ooo 

0. 

35 

o. 

1.000 

WITH  ADJACENT  DURATION 

AT 

ABOVE 

STATION 

32 

.527 

0. 

35 

.753 

.981 

COR'n'E 

L*TION 

COEFFICIENTS  OF  RECORDED  DATA 

FOR  30- 

DAY  DURATION 

STS 

32 

35 

WITH  SANE  DURATION 

TP 

1  .OCO 

0. 

35 

0. 

1  .000 

WITH  ADJACENT  duration 

AT 

ABOVE 

STATION 

32 

.690 

C. 

35 

0, 

,833 

RECORDED  ANO  RECONSTITUTED  DAT4 


STA 

YEAR 

PEAK 

1-DAY 

3-OAY 

10-DAY 

30-DAY 

32 

1905 

77100. 

71200. 

62000. 

51000. 

30830. 

32 

1906 

206000. 

185000, 

130000. 

83OC0. 

51000. 

32 

1907 

138000, 

133000. 

115000, 

653C0 . 

03670. 

3i 

1  908 

185000. 

167C0C. 

132000. 

65600. 

00130. 

32 

1909 

137000. 

122000. 

7  0  0  0  0  . 

b60C0. 

36130. 

32 

1930 

99o00 . 

95200. 

9C0C0, 

602C0, 

46100. 

STA 

YEAR 

PEAK 

l.CAY 

3-0  AY 

10-DAY 

30-DAY 

35 

1905 

2S357.E 

20007. E 

20362. E 

12398. E 

8837. E 

35 

i  9  06 

08000. 

32500. 

20300. 

1267c, 

7093. 

35 

1  007 

06000. 

32600. 

29270. 

16020. 

9570. 

35 

1948 

53000. 

00700. 

20870. 

129?  U . 

6890. 

»5 

1  949 

186.90. 

1 06  CO , 

10570. 

8090. 

5690. 

55 

1930 

23600. 

20100. 

15800. 

9800. 

6920. 
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FXHIRIT  4 


tREOUEN'CV  STATISTICS 

CF  RECORDED 

AND  1 

RECONSTITUTED 

DATA 

STA  ITEM 

PEAK 

1-DAY 

3-DAY 

10-DAY 

30-DAY 

32  MEAN 

5.123 

5.085 

4.985 

4.835 

4.621 

STD  DEV 

.161 

.154 

.142 

.083 

,076 

SKEW 

-.366 

-  .527 

-.375 

-  .266 

-1.088 

EQUIV  YKS 

6,0 

6.0 

6.0 

6.0 

6.0 

35  MEAN 

4.520 

4. "401 

4.257 

4.070 

3.873 

3T0  DEV 

.154 

.168 

.162 

.104 

.062 

SKEW 

-.176 

-.240 

-1.036 

.573 

-.012 

EDUIV  YRS 

5.4 

6.0 

5.8 

6.0 

5.8 

CORRELATION  COEFFICIENTS  OF  RECORDED  AND  RECONSTITUTED  DATA  FOR  PEAK  DURATION 
WITH  SAWE  curation 

WITH  ADJACENT  DURATION  AT  ABOVE  STATION 

:IENTS  OF  RECORDEO  AND  RECONSTITUTED  DATA  FOR  1-DAy  DURATION 
WITH  SAME  DURATION 

WITH  ADJACENT  DURATION  AT  ABOVE  STATION 


STA 

32 

35 

32 

1  .000 

.574 

35 

.57a 

1  .000 

32 

.955 

.475 

35 

.616 

.986 

?REI 

LATION 

COEFF 

3TA 

32 

35 

32 

1  .000 

.526 

35 

.526 

1.000 

32 

.955 

.616 

35 

.475 

.986 

CORRELATION  COEFFICIENTS  OF  RECOROEO  AND  RECONSTITUTED  DATA  FOR  3-DAY  DURATION 
OTA  32  35 

HITH  SAFE  DURATION 

32  l.OCO  .533 
35  .553  1.000 

WITH  ADJACENT  DURATION  AT  ABOVE  STATION 

32  .036  .608 

35  0.  .375 


CORRELATION  COEFFICIENTS  OF  RECORDEO  ANO  RECONSTITUTED  DATA  FOR  IO-OAy  DURATION 
STA  ??  55 

WITH  SANE  DURATION 

32  I. 000  0. 

35  0,  1.000 

HITH  ADJACENT  DURATION  AT  ABOVE  STATION 

32  .920  0. 

35  .365  .577 

CORRELATION  COEFFICIENTS  OF  RECORDED  ANO  RECONSTITUTED  DATA  FOR  30-DAY  DURATION 
STA  32  35 

WITH  SAME  CURATION 

32  1,000  0. 

35  0.  1.000 

WITH  ADJACENT  OURATICN  AT  ABOVE  STATION 

32  .7aa  o. 

35  0.  .815 


-4- 


EXHIRIT  « 


frequency  arrays 


STATION  32 


NO 

PLOT 

PEAK 

1-DAY 

3-OAY 

1C-DAY 

30-DAY 

1 

10. PI 

206000. 

*85000 . 

134000. 

85800. 

51000. 

2 

26. 55 

185000, 

16700C. 

122000. 

83400. 

461C0. 

3 

42.18 

130000. 

13300C. 

115000. 

66800. 

44130, 

4 

57.52 

137000. 

122000. 

90000 , 

65300. 

43670. 

5 

73. as 

99000. 

959CC. 

70400. 

64200. 

38130, 

6 

80.  OP 

77100. 

71200. 

62000. 

51000. 

30830. 

STATION 

35 

NO 

PLOT 

PEAK 

1-OAY 

3 -DAY 

10-DAY 

30-DAY 

1 

10.91 

53400. 

4030C. 

29270. 

16020. 

9570. 

2 

28 . 5? 

48400. 

32600. 

24670. 

12980. 

8837. E 

3 

*2.18 

46000. 

32500. 

24300. 

12870. 

7493. 

4 

57.82 

25357. E 

20407  .E 

20362. E 

12398. E 

6920. 

5 

73.45 

23600. 

20 1 00 . 

15800. 

9300. 

6890, 

6 

39.09 

18600, 

14600, 

10570. 

6C90. 

5690, 

AOOPTED 

sta 

frequency 

ITEM 

STATISTICS 

PEAK 

1-DAY 

3-DAY 

10-DAY 

30-DAY 

32 

HP  A  N 

5.123 

5.089 

4.984 

4.835 

4.621 

STO  DEV 

.159 

.153 

.133 

.106 

.066 

SKEW 

-.334 

-.366 

-.462 

-.599 

-.795 

INCRMT 

0. 

0. 

0. 

c. 

0. 

35 

MEAN 

4.520 

4.401 

4.297 

4.070 

3.873 

STO  DEV 

.193 

.172 

.153 

.113 

.078 

SKEW 

-.462 

-.437 

-.414 

-  .  365 

-.322 

INCRMT 

0. 

0. 

0. 

0. 

0. 

COMPUTED 

FREQUENCY 

CURVES 

STATION 

PLOT 

32 

EXP  PROB 

PEAK 

1-OAT 

3-DAY 

10-DAY 

30-DAY 

.01 

.74 

404931. 

350275. 

227309, 

127137. 

58610, 

.10 

1.84 

347604. 

304541. 

204648. 

119050. 

57022. 

1.00 

4.25 

284819, 

253204, 

177239. 

108179. 

54512. 

5.00 

9.65 

233536. 

210253. 

152569. 

97369, 

51625. 

10.00 

14.65 

209104. 

189451. 

140041  . 

.  91478. 

40894, 

30.00 

32.63 

162936. 

149501. 

114730. 

78768. 

43770. 

50.00 

50.00 

135505. 

125410. 

90620. 

70013. 

42b27. 

70.00 

67.37 

111592. 

10399c; . 

83702. 

614  44. 

39276. 

90.00 

85.35 

82252. 

77412. 

64306. 

4952b. 

34128, 

95.00 

90.35 

70S68. 

66693. 

56175. 

44236. 

31637. 

99.00 

95.75 

51767, 

49272. 

42524. 

34903. 

26879. 

99.90 

98.16 

35879. 

34363. 

30350. 

2bO 1 7 , 

21825, 

99,99 

99.26 

25794. 

24600, 

22260. 

1*1751. 

17860. 

STATION 

PLOT 

35 

EXP  PROB 

PEAK 

1-OAY 

3-DAY 

10-OAY 

30-UAY 

.01 

.79 

115212. 

77843. 

54915. 

25533. 

12892. 

.10 

1.94 

96922. 

67581. 

48195. 

23022. 

11952. 

1.00 

4.40 

80287 . 

55798. 

40443. 

20083. 

10832. 

5,00 

9. 84 

64601 . 

45790. 

33798. 

17503, 

9824. 

io.oo 

14.84 

57035. 

40914. 

30330. 

16205. 

9305. 

30.00 

32.72 

42704. 

31545. 

24167  , 

13601, 

8235. 

50,00 

50.00 

34203. 

2392ft, 

20200 , 

1 1  °44  . 

7529. 

70,00 

67.28 

27020, 

20980. 

16797. 

10399, 

6847. 

90.00 

85.16 

16« 2°. 

14956. 

12443. 

8  560, 

59u3  , 

95.00 

90.16 

15146. 

12582, 

10681 . 

7488. 

5480, 

9g,00 

95.60 

10111. 

8824. 

7816, 

5986. 

4716, 

99,90 

96.06 

6197. 

5753. 

5372, 

4586. 

3950. 

49,99 

99.2) 

3951, 

3889. 

3815. 

3604. 

3369. 
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EXHIRIT  4 


JULY  i 972  723-X6-L2350 
REGIONAL  FREQUENCY  COMPUTATION 
VERSION  DATE  -  AUGUST  21,  1979 


TEST  DAT* 

727-X6-L2350 

SAVE  STATIONS  FROM  PREVIOUS  JOB 


NDUR  I YRA 

ISKEH 

KEEP 

ICONV 

IPCHQ 

IPCHS  NSTAT 

NSMTH 

INCAO 

5  1945 

1 

2 

-0 

-0 

-0  -0 

-C 

-0 

REGIONAL  SKEW  COEFFICIENTS 

PEAK 

1-DAY 

3-DAY 

10-DAY 

30-DAY 

-.200 

-.400 

-.500 

-.600 

-.800 

STATION(S)  kept 

from  last 

RUN.  32. 

35, 

FREQUENCY  STATISTICS  OF  RECOROEO  DATA 

STA  ITEM 

PEAK 

1-DAY 

3-OAY 

10-DAY 

30-DAY 

33  MEAN 

3.941 

3.867 

3.781 

3.602 

3.443 

STD  DEV 

.158 

.137 

.148 

.090 

.083 

SKEW 

-.320 

-.599 

-.412 

-.371 

.180 

INCRMT 

0. 

0. 

0. 

0. 

0. 

YEARS 

5. 

5. 

s. 

5. 

5. 

FREQUENCY  STATISTICS  AFTER  ADJUSTMENT  WITH  A  LONG  TERM  STATION 


STA  ITEM 

PEAK 

1-OAY 

3-OAy 

10-DAY 

32  MEAN 

5.123 

5.089 

4.984 

4.835 

STD 

DEV 

.161 

.154 

.1«2 

.083 

SKEW 

-.200 

-.400 

-.500 

-.600 

INCRMT 

0. 

0. 

0. 

0. 

EQUIV 

yrs 

6.0 

6.0 

6.0 

6.0 

35  MEAN 

4.520 

4.401 

4.297 

4.070 

STD 

OEV 

.194 

.168 

.162 

.104 

SKEW 

-.200 

-.400 

-.500 

-  .600 

INCRMT 

0. 

0. 

0. 

0. 

EQUIV 

YRS 

6.0 

6.0 

6.0 

6.0 

33  l 

MEAN 

3.921 

3,854 

3.776 

3.593 

STD 

OEV 

.155 

.131 

.135 

.088 

SKEW 

-.200 

-.400 

-.500 

-.600 

INCRMT 

0. 

0. 

0. 

0. 

EQUIV 

YRS 

5.6 

5.7 

5.8 

5.5 

CORRELATION 

COEFFICIENTS  OF 

RECOROEO  DATA 

FOR  PEAK  DURATION 

STA  32 

35 

33 

WITH  SAME 

DURA'.  ION 

32  1,000 

.5R9 

.828 

35  .599 

1  .000 

.824 

33  .828 

.824 

1.000 

WITH  ADJACENT  DURATION 

AT  ABOVE 

STATION 

32  .996 

.520 

.701 

35  .652 

.985 

.911 

33  .873 

.730 

.955 

CORRELATION 

coefficients  of 

RECOROEO  DATA 

FOR  1-1 

DAY  DURATION 

STa  32 

35 

33 

WITH  SAME 

DURATION 

32  1.000 

.579 

.768 

35  .579 

1.000 

.845 

33  .768 

.845 

1.000 

WITH  ADJACENT  DURATION 

AT  ABOVE 

STATION 

32  .998 

.652 

.873 

35  .520 

.’85 

.730 

33  .701 

.’ll 

.955 

30-DAY 

4.621 

.076 

-.000 

0. 

6.0 

3.673 

.082 

-.800 

0. 

6.0 


3.449 

.080 

-.800 

0. 

5.4 
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EXHIBIT  4 


CORRELATION  COEFFICIENTS  QF  RECORDED  0*T*  FOR  3-OAT  DURATION 

STA  32  35  33 

WITH  SAME  DURATION 

32  J.000  .588  .876 

35  .588  1.000  .729 

33  .876  .728  1.000 

WITH  ADJACENT  DURATION  AT  ABOVE  STATION 

32  .050  .867  .963 

35  0.  .857  ,679 

33  .623  .781  .979 


CORRELATION  COEFFICIENTS  OF  RECORDED  OATA  FOR  10-0AY  DURATION 

STA  32  35  33 

HITH  3 AME  DURATION 

32  1.000  0.  ,297 

35  0.  1.000  ,708 

33  .297  .708  1,000 

WITH  AOJlCENT  DURATION  AT  ABOVE  STATIC;; 

32  .8 28  0.  .335 

35  .383  .968  .666 

33  .050  ,783  .973 


CORRELATION  COEFFICIENTS  OF  RECORDED  DATA  FOR  30-DAY  DURATION 

STA  32  35  33 

WITH  SAME  DURATION 

32  1.000  0.  .586 

35  0.  1.000  .183 

33  .586  .183  1.000 

WITH  ADJACENT  DURATION  AT  ABOVE  STATION 

32  .690  0.  .769 

35  C.  .805  0. 

33  0.  .706  .920 


RECORDED  AND  RECONSTITUTED  DATA 


STA 

year 

peak 

i-dat 

3-DAY 

10-DAY 

30-DAY 

33 

1995 

5530. 

5040, 

4100. 

3320. 

2270. 

33 

1996 

13300. 

9560. 

7700. 

4840. 

3150. 

33 

1997 

10300. 

9360. 

8530. 

4850. 

3540. 

33 

1998 

1 0300  , 

0840. 

6930. 

4230. 

2790, 

33 

1949 

6470. 

5900. 

4300. 

3120. 

2330, 

33 

;9so 

6669. E 

6246. E 

6157. E 

«1S1.E 

3131, E 

FREQUENCY  STATISTICS 

OF  RECORCEO 

AND 

RECONSTITUTED 

DATA 

STA  ITEM 

PE»K 

1-0  AY 

3-DAy 

10-DAY 

30-DAY 

33  MEAN 

3.921 

3.855 

3.782 

3.605 

3.452 

3T0  OE V 

.149 

.126 

.132 

.081 

.077 

SKEW 

.173 

-.185 

-.462 

-.513 

-.256 

EQUIV  YRS 

5.9 

5.9 

5.9 

5.9 

5.9 

CORRELATION  COEFFICIENTS  QF  RECORDED  And  RECONSTITUTED  DATA  FOR  PEAK  CURATION 

STA  32  35  33 

WITH  SAME  DURATION 

32  l.POO  .579  .855 

35  .579  1.000  .853 

33  .855  .853  1,000 

WITH  ADJACENT  DURATION  AT  ABOVE  STATION 

32  .995  .475  .739 

35  .616  .986  .907 

33  .087  .777  .95’ 
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EXHIBIT  4 


CORRELATION  COEFFICIENTS  OF  RECORDED  AND  RECONSTITUTED  DATA  FOR  1-DAY  DURATION 


STA 

32 

35 

33 

HITM 

SAME  DURATION 

32 

1  .004 

.526 

.795 

35 

.326 

1.000 

.864 

33 

.795 

.864 

1  .000 
WITH 

ADJACENT  DURATION  AT  ABOVE  STATION 

32 

.995 

.616 

.887 

35 

.975 

.936 

.777 

33 

.739 

.907 

.959 

correl 

.  A  T  I  0  N 

COEFFICIENTS 

OF  RECORDED  AND 

RECONSTITUTED  DATA 

FQR 

STA 

32 

35 

33 

WITH 

SAME  DURATION 

32 

1.000 

.558 

.904 

35 

.553 

1.000 

.647 

-j  2j 

.904 

.647 

1.000 

WITH 

adjacent  duration 

AT  ABOVE  STATION 

32 

.836 

,848 

.964 

35 

0. 

.875 

.686 

33 

.604 

.761 

,9ill 

CCRREL 

.  A  T I  ON 

COEFFICIENTS 

OF  RECORDED  AND 

RECONSTITUTED  DATA 

FOR 

STA 

32 

35 

33 

WITH 

same  duration 

32 

1.000 

0. 

.375 

35 

0. 

1.000 

.615 

33 

.375 

.615 

1  .000 

WITH 

ADJACENT  DURATION 

AT  ABOVE  STATION 

32 

.820 

C. 

.498 

35 

.3b3 

.977 

.561 

33 

."61 

.709 

.973 

- 

•  - 

CORRELATION 

COEFFICIENTS 

OF  RECORDED  ANQ 

RECONSTITUTED  DATA 

FOR 

STA 

32 

35 

33 

KITH 

SAME  CURATION 

32 

1 .000 

c. 

.740 

35 

0. 

1.000 

0. 

33 

.740 

0. 

1.000 

WITH 

ADJACENT  DURATION 

AT  ABOVE  STATION 

32 

.744 

0. 

.751 

35 

0. 

.819 

0. 

33 

0. 

.365 

.923 

3-DAY  DURATION 


10-DAY  DURATION 


30-day  duration 


FPE3UENCY 

ARRAYS 

STATION 

33 

NO  FLCT 

peak 

l-DAY 

1  10.91 

13300. 

956C. 

2  26.55 

10300. 

9360. 

3  42.18 

10300. 

884C. 

4  57.82 

6669. E 

6246. E 

3  73.35 

6470. 

540C. 

6  89,09 

5530. 

5040. 

3-OAY 

1C-DAY 

30-DAY 

8530. 

4850. 

3540. 

7700. 

4840. 

3150. 

6930. 

4230. 

3131. E 

6157. E 

4151.  E 

2790, 

4300, 

3320, 

2330. 

4100. 

3120. 

2270. 

-8- 


FXHI8IT  4 


ADOPTED 

sta 

FREQUENCY 

ITEM 

statistics 

PEAK 

1-DAY 

3-DAY 

10-DAY 

30-DAY 

33 

MEAN 

3.921 

3.855 

3.782 

3.605 

3.952 

STD  OEV 

.195 

.139 

.123 

.099 

.070 

SKEW 

-.200 

-,«C0 

-.500 

-.600 

-.800 

INCRMT 

0. 

0. 

0. 

0. 

0. 

COMPUTED  FREQUENCY  CURVES 
STATION  33 


PLOT 

EXP  PROS 

peak 

l-OAY 

3-DAY 

10-DAY 

30-DA 

.01 

.76 

25133. 

17617. 

13090. 

7011. 

9096 

.10 

1.88 

21266. 

15677. 

11953. 

6611. 

3932 

1.00 

9.30 

17255. 

13919. 

10535. 

6070. 

3750 

5.00 

9.72 

19127. 

11995. 

9218. 

5525. 

3591 

10.00 

19.72 

12677. 

10963. 

8539. 

5226. 

3916 

30.00 

32.66 

9995. 

8519. 

7122, 

9571. 

3118 

50.00 

50.00 

8939. 

7306. 

6201. 

9115. 

2892 

70.00 

67.39 

7079, 

6196. 

5331. 

3662. 

2652 

90.00 

85.28 

5910, 

9770. 

9173. 

3020. 

2285 

95.00 

90.28 

9793  . 

9177. 

3678. 

2730. 

2108 

99.00 

95.70 

3659. 

3185. 

2632. 

2209. 

1773 

99,90 

98.12 

2706. 

2302. 

2059. 

1699. 

1921 

99.99 

99.29 

2078. 

1713. 

1539. 

1326. 

1198 
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EXHIBIT  4 


JULY  197  2  723-X6-L2350 
REGIONAL  '■requency  COMPUTATION 
VERSION  DATE  -  AUGUST  21,  1979 


TEST  DATA 
723-X6-L2350 
STATISTICS  FURNISHED 


NDUR 

IYRA 

I  SKEW 

KEEP 

S 

1945 

1 

-0 

ADOPTED 

FREQUENCY 

STATISTICS 

STA 

ITEM 

PEAK 

1-DAY 

32 

MEAN 

5.123 

5.089 

STO  OEV 

.159 

.153 

SKEW 

-.200 

-.400 

INCRMT 

0. 

0. 

35 

MEAN 

4.518 

4.408 

STO  DEV 

.196 

.177 

SKEW 

-.200 

-.400 

INCRMT 

0. 

0. 

REGIONAL  SKEW  coefficients 

peak 

-.200 

INPUT  FREQUENCY  STATISTICS 

1-DAY 

-.400 

STA  ITEM 

PEAK 

1-DAY 

32  MEAN 

5.123 

5.089 

STO  OEV 

.159 

.153 

SKEW 

-.334 

- .  366 

INCRMT 

0. 

0. 

EQUIV  VAS 

6.0 

6.0 

35  MEAN 

4.518 

4.008 

STO  OEV 

.196 

.177 

SKEW 

-.278 

-.168 

INCRMT 

0. 

0. 

Equiv  VRS 

5.6 

6.0 

iconv 

IPCHQ 

IPCHS 

NSTAT 

-0 

-0 

•0 

2 

3-OAY 

10-DAY 

30-DAY 

4.984 

4.835 

4.621 

.133 

.106 

.066 

-.500 

-.600 

-.800 

0. 

C. 

0. 

4,267 

4.052 

3.843 

.153 

.117 

.082 

-.500 

-.600 

-,eoo 

0. 

C. 

0. 

3-DAY 

1C-DAY 

30-DAY 

-.500 

-.600 

-.600 

3-CAY 

1C-DAY 

30-DAY 

4.984 

4.635 

4.621 

.133 

.106 

.066 

-.462 

-.599 

-.795 

0. 

c. 

0. 

6.0 

6  .0 

6.0 

4.267 

4.C52 

3.843 

.153 

.117 

.082 

-.027 

.168 

.396 

0. 

0. 

0. 

6.0 

5.8 

5.9 

NSMTH 

-1 


INCAD 

-0 
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EXHIBIT  A 


computed  frequency  curves 


station 

PLOT 

.01 

.10 

1.00 

5.00 

10.00 

30.00 

50.00 

70.00 

90.00 

05.00 

05,00 

95.90 

99,99 


32 

EX?  PRD8 
.74 
1.84 
4.25 
9.65 
14.65 
32.63 
50.00 
67.37 
8  5.35 
90.3S 
55.75 
*6.16 
99. S6 


PEAK 
4461 15, 
371280. 
295034. 
236862. 
210267. 
161945. 
134367. 
110754. 
32495. 
71366. 
53600. 
38329. 
28829. 


1-DAY 

342977, 

300235. 

251317. 

209670. 

169276. 

149725. 

125646. 

104117. 

77257. 

66397. 

48734. 

33646. 

24015, 


3-DAV 

222370. 

201503. 

175691. 

152007. 

139799. 

114666. 

96857. 

93904. 

64335. 

56098. 

02243. 

29892. 

21711, 


station 

35 

PLOT 

EXP  PR08 

.01 

.7  r 

.10 

i.’i 

1.00 

4.34 

5.00 

9.77 

10.00 

14.77 

30,00 

32.69 

50.00 

50.00 

70.00 

67.31 

90.00 

85.23 

95.00 

<90.23 

99,00 

95.66 

99.90 

96.09 

99.99 

99.23 

PEAK  1-DAY 

146877.  83998. 

117125.  72011. 

80243.  58619. 

67301.  47535. 

581  13,  42228 . 

02118.  32196. 

33460.  26287. 

?63&7 .  21150. 

10338.  1497fa. 

15344,  12569, 

10777,  8788. 

7174.  5725 . 

5018.  3876. 


3-DAY 

48381. 

03196. 

36694, 

31233. 

28365. 

22633. 

19040. 

15766. 

11616. 

9922. 

7160. 

4809. 

3329. 


10-DAY 

30-DAY 

127450. 

56504. 

119328. 

56939. 

108407, 

54454. 

97546. 

51588. 

91627. 

49866. 

78837, 

45759, 

70064. 

42623. 

61459. 

39277. 

49497. 

34132, 

44191. 

31641, 

34633. 

26879. 

25933. 

21818. 

19663. 

17e48. 

10-DAY 

30-UAy 

22406. 

10563. 

20637. 

10233. 

18742. 

9661  . 

16661 . 

9052. 

15567. 

8678. 

13167. 

7799. 

11577. 

7141. 

10016. 

6451. 

7839. 

5418. 

6961. 

4931  . 

5353. 

4027. 

3665. 

3108, 

2648. 

2421  . 

EXHIBIT  4 


EXHIBIT  5 


DEFINITIONS  723-X6-L7350 


AA(I) 

AB(I) 

ABS 

ALOG 

ANYR(l,K) 
ANYRS 
AV ( I j  K) 

AVGSK 

b(k) 

BB 

BC 

BLANK 

CB 

CC 

CROUT 

DQ(I,K) 

DTRMC 

E 

I 

IA 

ICORL 

ICSE 


II 

INCAD 

INDC 

IPCHQ 

IPCHS 

IPREV 

IRCRD(j) 

IRATO 

ISKEW 


-  First  half  of  description  for  duration  I 

-  Second  half 

-  Computer  library  function  for  absolute  value  of 
number 

-  Computer  library  function  for  natural  logarithm 

-  Number  of  years  of  data  for  station  K  and  duration  I 

-  Number  of  years  of  data  in  study 

-  Mean  logarithm  (or  sum  of  logarithms)  for  station  K 
and  duration  I 

-  Average  regional  skew  coefficient 

-  Regression  coefficient  for  variable  (k) 

-  Regression  coefficient 

-  Regression  coefficient 

-  Symbol  to  identify  recorded  data 

-  Regression  constant 

-  Regression  constant 

-  Program  subroutine  to  solve  simultaneous  equations 

-  Increment  added  to  all  flows  for  duration  I  at  station 
K  to  preclude  infinite  negative  logarithms 

-  Multiple  determination  coefficient 

-  Symbol  to  identify  reconstituted  data 

-  Index  for  duration 

-  Indicator  in  column  1  of  first  card  for  each  job 

-  Indicator,  when  positive  calls  for  computation  of 
correlation  coefficients 

-  Indicator,  case  number  specifing  cause  for  no  independent 
variables  in  estimation  equation 

+1  indicates  no  flows  found  for  correlation 
+2  indicates  all  correlations  were  zero 

-  Index  associated  with  I 

-  Indicator,  positive  value  calls  for  adjustment  of 
increment  to  reduce  skew  coefficient 

-  Indicator  positive  when  correlation  coefficient  has 
been  changed 

-  Indicator,  when  positive  calls  for  punching  recorded 
and  reconstituted  flows 

-  Indicator,  when  positive  calls  for  punching  statistics 

-  Order  number  in  regression  equation  of  adjacent  duration 

-  Indicator  blank  when  no  record  at  all  stations  in  year  J 

-  Indicator,  when  positive  calls  for  reading  conversion 
ratios 

-  Indicator  when  positive  calls  for  reading  skew  coefficients 


EXHIBIT  5 


ISTA(K) 

ISTAN 

ISTN 

ISTY 

ITEMP 

ITMP 

ITP 

IX 

IXX 

IYR 

IYRA 

J 

JA 

JX 

K 

KDUR 

KEEP 

KEPT(K) 

KRCRD 

KSTA 

KX 

KYRS 

L 

LA 

LTRA 

LX 

M 

MM 

N 

NCAB(I,K) 

NDUR 

NINDP 

NLOG(l,K) 

NSMTH 

NSTA 
NS  TAT 

NSTAX 

NSTXX 

NVAR 

NYDIF 

NYRS 


Identification  number  for  station  K 
Station  number 

Array  of  station  sequence  by  length  of  record;  longest 
record  first 

Array  of  station  record  lengths  used  to  build  ISTN  array 

Temporary  variable 

Temporary  variable 

Temporary  variable 

Index  associated  with  I 

Argument  for  random  number  function 

Year  number 

Number  of  earliest  year  of  record 
Year  index 

Index  associated  with  J 
Index  associated  with  J 
Station  index 

Dimension  limit  for  durations 

Number  stations  to  keep  from  immediately  previous  job 
Station  numbers  kept  from  immediately  previous  job 
Indicator,  when  positive  a  complete  record  exists  for 
all  stations 

Dimension  limit  for  stations 
Index  associated  with  K 
Dimension  limit  for  years 
Subordinate  station  index 
Index  associated  with  L 
Letter  A  for  testing  IA 
Index  associated  with  L 
Sequence  index 
Index  associated  with  M 
Temporary  counter 

Number  of  cross  products  for  station  K  and  duration  I 
Number  of  durations  in  study 

Number  of  independent  variables  in  correlation 
Number  of  values  for  station  K  and  duration  I 
Indicator,  zero  or  positive  value  causes  smoothing  of 
statistics 

Number  of  stations  in  study 

Number  of  stations  for  which  statistics  (instead  of 
flows)  are  supplied 
Twice  NSTA 

Number  of  stations  kept  from  previous  job  incremented  by  1 
Total  number  of  variables  in  correlation 
Indicator,  when  positive  a  difference  in  record  length 
exists  between  new  data  and  data  from  previous  job 
Number  of  years  in  study 
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EXHIBIT  5 


P(I) 

PLTT(J) 

Q(M,K) 

QM(I) 

QR(M,K)  - 

R(K,K+l)  - 

RA(I,K,L)  - 

RMAX 

RMIN 

SA 

SAA 

SAB 

SAC 

SB 

SC 

SD(I,K) 

SDA 

SDB 

SIN 

SKEW(I,K)  - 

SKW(l) 
SQA(I,K)  - 

sqb(i»k)  - 
suma(i,k)  - 

SUMB(I,K)  - 
T 

TEMP 

TMP 

TMPA 

TMPB 

TMPP 

TP 

X(K) 

XINCR(I,K)- 
XPAB(I,K)  - 

xq(!) 


Exceedence  frequency  coordinate  or  ratio  to  convert 

flows  to  average  rates 

Plotting  position  for  event  number  J 

Flow  or  logarithm  for  station  K  and  sequence  number  M 

Flow  for  current  station  and  year  and  for  duration  I 

Indicator  whether  Q(M,K)  is  recorded  or  reconstituted 

Covariance  array  for  multiple  regression  equation 

Correlation  between  stations  K  and  L  for  duration  I 

Maximum  consistent  correlation  coefficient 

Minimum  consistent  correlation  coefficient 

Sum  of  mean  logarithms  for  various  durations 

Sum  of  squares  of  mean  logarithms 

Sum  of  cross  products  of  mean  logarithm  and  standard 
deviation 

Sum  of  cross  products  of  mean  logarithm  and  skew 
coefficient 

Sum  of  standard  deviations  for  various  durations 
Sum  of  skew  coefficients  for  various  durations 
Standard  deviation  (or  sum  of  squares)  for  station  K 
and  duration  I 

Standard  deviation  of  short  record  station 
Standard  deviation  of  long  record  station 
Computer  library  function  for  sine 
Skew  coefficient  (or  sum  of  cubes)  for  station  K  and 
duration  I 

Specified  skew  coefficient  for  duration  I  at  all  stations 
Sum  of  squares  of  logarithms  in  correlation  for  station 
K  and  duration  I 

Sum  of  squares  of  logarithms  at  related  station  in 
correlation  with  station  K  for  duration  I 
Sum  of  logarithms  in  correlation  for  station  K  and 
duration  I 

Sum  of  logarithms  at  related  station  in  correlation 

with  station  K  for  duration  I 

Large  number  denoting  missing  record 

Temporary  variable 

Temporary  variable 

Temporary  variable 

Temporary  variable 

Temporary  variable 

Temporary  variable 

Independent  variable  related  to  station  K 
Increment  for  DQ  in  skew  coefficient  adjustment  routine 
Sum  of  cross  products  of  logarithms  for  station  K  with 
related  station  for  duration  I 
Temporary  flow  array 
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'>  n  o  o  o  n  n  o 


SOURCE  PROGRAM 


723-X6-L7350  REGIONAL  FREQUENCY  COMPUTATION,  HEC,  JULY  1973  1001 

LIBRARY  SUBROUTINES  USED— ALCG,slN,A0s  1002 

PROGRAM  SUBROUTINES  cocut, RnGKn— see  COMMENTS  In  RNGEn  1003 

REFERENCE  TO  tape  ?  AT  960+1,1170+0  1004 

indexes  i=duf!At:on  j-year  k=station  l=related  sta  h=sequence  no  1005 

1006 

DIMENSION  1007 

m(8),AB(8),ANYR(8,10),AV(8,10)fB(10),DQ(8,iO)#  lOoti 

2!RCRl-(100),  1ST  A  (10),  ISTM10),  1ST  Y  (10)  ,  KEPT  (10)  ,NCA8  (8,10,20),  1009 

3NL0S;E, 10) ,P(S) ,PLTT( IOC) ,Q (000,10) ,0M  (400) ,QMIN (8,10) ,  10 10 

03*?(a00.10),R(lC,ll),RA(e,10,20),SD(8,lC),SKEH(8,10),SKW(8),  1011 

SSQA (3, 10,20) » S  0  B (fl, 10,20 ,3UMA(8, 10,20) ,SUMO(e,  1C,20) ,X(flOO) ,  1 0 1 2 

6XINCR(8,10),XPAD(li,10,2C),XQ(B)  l  0  1  3 

COMMON  0 TRMC , N I NDP , 3  l0i« 

DATA  LTRA/1HA/.BLANK/1H  /,E/1HE/  1015 

KSTA=10  1016 

KDUR=8  1017 

XYSS=50  1018 

1C  FORMAT (IX,  17,918)  1019 

20  FORMAT (1X,F7.0»9F8.0)  1020 

30  FORMAT (Al,A3»9A4,10A4)  1021 

40  FORMAT (IX, A3, 9Aa,lOAA)  1022 

50  FORMATHHi)  1023 

60  FORMATflX, 17,18, 8F8.0)  1029 

70  FORMAT  (2X,A5,A4,F9.3)  1025 

80  FORMAT  (lX,2Aa,F9.3)  1026 

00  90  K=1,*$TA  1027 

90  1STA(K)=-1  1026 

IYRSV=0  1029 

WASTE  CAROS  UNTIL  AN  A  IN  COL  1,  FIRST  TITLE  CARO  1030 

**  CARD  A"1  **  1031 

100  READ(5,3C)IA,  (QS(J, l),Jsl,20)  1032 

ir  (IA.NE.LTRA)  GO  TO  ICO  1033 

C  **  CARO  A-2,3  *•  1034 

OCAD(5,aO) ((QR(J,K),J=1,20),K=2,3)  1035 

0  **  CARD  B  **  1036 

RE A 0(5, 10)  N0UR,IYRA, ISKEW,KEEP, ICONV, IpCHQ, IPCHS, NSTAT, NSMTH,  1037 

1INCAD  1038 

c  TERMINATE  WITH  4  BLANK  CARDS,  AN  A  IN  CCL  1  OF  FIRST  1039 

IF(NCUR.LE.Q)  STOP  1040 

HR  I  T£  (6 , 30)  1041 

WRITE(6, 1 10)  10^2 

no  F0RmAT(1X,3O(1H*)/10H  JULY  1972,9x,12H723-xfc-L2350/9H  REGIONAL,  1043 

s  22w  frequency  computaticn/3ih  version  date  -  august  21,  1979/  104a 

s  ix,30(tM+)///)  lOaai 

WRITE (6, 40)  ((QR(J, K),J= 1,20), Ka 1,3)  1045 

IF (NDUR ,LE .KOURlGO  TO  1«0  1046 

J  20  WRITE (6, 130)  NSTA,NOUR,NYRS  1047 

130  fmrmaT(/19H  DIMENSION  EXCEEDED  , 5 X , 5HNST A  = , 1 3, 5X , 5HN0uR= , 12, 5X , 5HN  1048 
lVRSs.la)  1049 

GO  TO  100  1050 

140  KR  T  TE ( 6 , 150)  NDIR,IYSA,ISKEW, KEEP, ICONV, IPCHQ, IPCHS, NSTAT, N3HTH,  1051 
11'lCAD  1052 

150  FoRMAT(/GX,4HnOur,6x,4HIyra,5x,5HI3KEw,6X,4HKEEP,5X,5HICONV,SX,5HI  1053 
1PCH3,5X,5HIPCHS,5X,5HNSTAT,5X,5HN5MTH,5X,5HINCA0,/10H0)  1054 

C  **  CARO  C  **  1055 

READ (5,40) ( A  A ( I ) »  A  3  ( I ) ,1=1, NO UR)  1056 

IF(ISKEw.LF.O)30  TO  200  1057 

A'/GSKsO,  1058 

C  **  CARD  D  •*  1059 

REA0(5,20) (SKW(I),I  =  1,N0UR)  _  1060 

WRITE (6 » 1 60)  1061 

ISO  F0RMATC/27H  REGIONAL  SKEW  COEFFICIENTS)  1062 

WRITE  (6, 170)  (A4(I),A8(I),T  =  1,NDUR)  1063 

17C  F[iRMATC20X.AJ,A4,7(3X,2A4))  J064 

W°ITE  (6, 180)  (3KW(I),Isi,N0UR)  1065 

ISO  FOR'MAT(16X,IOFU.3)  1046 

DO  190  1=1, NOUS  1067 

1 90  AvGSK=AVGSK+SKW(I)  1068 
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TtMP=NDU«  1069 

Avr,SK  =  AVGSK/TEHP  1070 

200  T  =  c!<?°:):;<)99.  107  i 

I:<Xs0  1072 

X7RA”IYR*«*1  1073 

N3TA=NSTAT  1074 

XF(NSTAT.GT.lO)  NSTAslO  1075 

NSTXX=1  1076 

IFCNSTAT.GT.O)  GO  TO  300  1077 

NSTA=0  1076 

INDC=0  1079 

NYOIF-0  1060 

C  INITIATE  -1,  NC  RECORO  FOR  ALL  FLOWS  1061 

ITd=KDUR*KYRS  1082 

DO  210  K=1,KSTA  1063 

DO  210  NSI.ITP  1084 

0R(N,K)=(-1.)  1065 

210  CONTINUE  1066 

IF(KEEP.LE.O)  go  TO  300  1067 

C  SAVE  STATIONS  FROM  PREVIOUS  RUN  IF  NECESSARY  1068 

C  **  CARD  E  **  1089 

READ  (5,10)  (KEPTCI), 1  =  1, KEEP)  1090 

KRITE(6,220.  (KEPTCI), IsirKEEP)  1091 

220  F0RMATC/31H  STATION(S)  KEPT  FROM  LAST  RUN  , 14 ( 1 H , 16) /31X6 ( 1H, 16) )  1092 
00  ?F0  K=1,KSTA  1093 

DO  270  Ls  1 , KEEP  1094 

IF  (KEPT(L)  . NE.ISTACK)  ''  GO  TO  270  1095 

INOC-I  1096 

NSTA=NSTAM  1097 

!STA(NSTA)=ISTA(K)  1098 

DO  230  Is!, NDUR  1099 

NLOG(I,NSTA)=0  1100 

0(3(1, NSTA)anQ(I,K)  1101 

XINCR(I,NSTA)SXINCR(I,K)  1102 

230  CONTINUE  1103 

M=0  1104 

ITM?=IYRSV-IYRA  1105 

►TI=ITMP*NDUR  1106 

ITP=IYRA-IyR3V-M  1107 

Ir(ITP.lE.O)  ITP=1  1108 

IF(MM.GE.O)  GO  TO  240  1109 

Ms-MM  1110 

MM=0  llil 

240  DO  260  J=ITP,NYRS  1112 

00  250  I=1,NCUR  1113 

M  =  *M  1114 

koshm+i  1115 

IF(IRCRD(J).LE.O)  GO  TO  250  1116 

TMP=0(M,K)  1117 

IF(TMP.GE.T)  GO  TO  250  1118 

QR(MM,NSTA)STMP  Ul9 

Nl_0G(I,NSTA)3NL0G{I*N3TA)*l  1120 

250  CONTINUE  1121 

260  CONTINUE  1122 

GO  TO  280  1123 

270  CONTINUE  1124 

280  CONTINUE  1125 

IF(ITMP.NE.O)  NTOIFal  1126 

NYR3=NY9S* ITMP  1127 

NSTXX=NSTA^1  11271 

IFCN3T4.Eo.KEEP)  GO  TO  300  1128 

IT°sKEEP-NSTA  1129 

WRITE(6,24C)  ITP  1130 

290  F OPn A  1 ( i 7 H  onT  ABLE  TO  F IND, 13 , 9HST ATIONS  )  1131 

KEcF=NSTA  1132 

300  IF ClvDC.LT.l)  NYR3=0  1133 

TF(ICONV.LE.O)  GO  TO  320  1134 

C  m  CARD  F  **  1135 

RE»D  (5,20)  (P(I), 1*1, ND'jR)  1136 

WRITE  (6,310)  1137 
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310 

FORMAT  C/30H  PATIOS  TO  OBTAIN  RATE  OF  FLOW) 

me 

WRITE (6, 170)  (AACIT  ,  A8(I)  ,I=1,NDUR) 

1139 

KSlTE(6»180) (P(I),I=t,NDUR) 

1140 

GO  TO  340 

1141 

320 

OC  330  I=I,NOUR 

1142 

PCI)=  1. 

1143 

320 

continue 

1144 

c 

SCT  CONSTANTS 

1145 

340 

IF(nSTAT.GT.O)  GO  TO  2140 

1146 

00  350  K=NSTXX,KSTA 

1147 

ISTA(K)=-1 

114B 

350 

continue 

1U9 

IYRSV=IYRA 

1150 

itp=kdur*kyrs/ndur 

1151 

DO  390  K=1,KSTA 

1152 

DO  330  1=1, HOUR 

1153 

IF(K.LT.NSTXX)  GO  TO  360 

1154 

NLGG(I,K)=0 

1155 

DO(I,K)=0. 

1156 

360 

DO  370  J=1 , ITP 

1157 

N=NDUR*(J-1)*I 

1158 

Q(N, K)=GR(N,K) 

1159 

370 

continue 

1160 

380 

CONTINUE 

1161 

390 

CONTINUE 

1162 

c 

* 

*  *  *  *  •  READ  and  process  ONE  STATION-YEAR 

DF  DATA 

*  *  *  1163 

c 

*•  CARD 

G  *• 

1164 

400 

read  (5,6  0) 1ST an, I YR, (QM(I) , I  =  1,NDUR) 

1165 

c 

**  CARD 

H  ** 

1166 

c 

BLANK  CARD  INDICATES  END  OF  FLOW  DATA 

1167 

IF(ISTAN.LT.1)G0  TO  470 

1168 

l'F(NSTA,LT.DGO  TO  420 

1169 

DO  410  K  =  1  , NST A 

1170 

c 

IOENTIFY  STATION  SUBSCRIPT 

1171 

IF(ISTAN.EQ.ISTA(K))GO  TO  430 

1172 

CIO 

CONTINUE 

1173 

420 

NSTA  =  NSTA-*1 

1174 

c 

ASSIGN  SUBSCRIPT  TO  NEW  STATION 

1175 

IF(NSTA.GT.KSTA)  GO  TO  120 

1176 

K=NSTA 

1177 

ISTA(K)=I3TAN 

1178 

c 

ASSIGN  SUBSCRIPT  TO  YEAR 

1179 

430 

J=IYR-IYRA 

1180 

IF(NYR5.LT. J)NYRS=J 

1181 

IF  (  j.GT .0) GO  TO  450 

1182 

WRITE(6,4uo)IYB 

1183 

440 

FORMAT ( / 1 8H  UNACCEPTABLE  YEAR  15) 

1184 

GO  TO  100 

1185 

c 

STORE  FLOWS  IN  STATION  AND  DURATION  ARRAY 

1186 

450 

M=(J-l)*NnuR 

1187 

DO  46C  I=1,N0UR 

1188 

M  =  M*1 

1169 

IF(GM(I).LE.C-1.))  GO  TO  460 

1190 

NLOG(I,K)=NLOG(I,K)*1 

_ 

1191 

D0(I,K)=D0(I,K)*QM(I) 

1192 

0(M,K)=GM(I) 

1193 

460 

continue 

1194 

GO  TO  400 

1195 

470 

IF(NYRS»N0UR.GT.KYRS*KCL'R)  GO  to  120 

1196 

c 

* 

»  *  *  *  *  COMPUTE  frequency  statistics  *  * 

ft  ft  ft  ft 

•  *  *  1197 

WRITE  (6 , 4f)0) 

1198 

480 

F0RMATC/38H  FREQUENCY  STATISTICS  OF  RECORDED  DATA  ) 

1199 

WRTTE(6,490) (At (I) ,A8(I) ,I=l,NDuR) 

1200 

490 

F()RMAT(SX,t2HSTA  ITEP  3x,A3,  A4,7  (3x,2A4)) 

1201 

00  500  J  =  1 »  NYRS 

1202 

500 

I RCRO ( J ) =0 

1203 

K  RC  R0  =  1 

1204 

ICDRL=1 

1205 

IF(NDUR.EQ.I.AN0.N3TA.EG.I)  ICORL=0 

1206 
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InDC=0 

oo  no  k*i,nsta 

t«pp=t 

xm!n=T 

OO  530  I=i,N9UR 
N  =  0 

IFCK.LT.NSTXX)  GO  TO  550 

TE^P'T 

M=I-NCUR 

O'  510  J=1,NYRS 

Msll  +  NOL'R 

TNf'  =  3(H,K) 

IF  (TtfP.LE  . (-1  ») )  GO  TO  510 
IFCTMP.LT.TEMP)  TEMPSTMP 
510  CONTINUE 

G!‘lN(I,K)=TEMP 

IF  (TEf'P.LT.THPP)  TMPPSTEMP 

TtHO=NLCG(I,K) 

T F  (TEmP.LT.0,1)  GO  TO  520 
Do(I,K)=0Q(I,K)*.001/7EP? 

IF fOO(I,K) ,L7.  .001)  DO  (  I ,  K )  =  .001 
TEMP=(GHIN(I,K)*DrJ(I.K))/OQ(I,K) 
Ic(T£MP.i-T.XMIN)  XKINsTEUP 
520  CONTINUE 

Dll  5 '10  1  =  1, NOUS 
IF(SUCG(I,X).LE.O)  GO  TC  530 

XINCR(I»K)=XHlN/16.*C0(I,K) 
IF(XINCH(I,K).LT..01)  XINCR(I,K)=,01 
530  IF(TMPP.GT,O..AND.lNCAO.LE.O)  OQ(I»K)eO, 

540  CONTINUE 
550  00  550  I=1,N0UR 
ANYR(I,K)=0. 

AvCIf K)=0. 

30(1, K)=0. 

S^F.  k'I»K)=0, 

56o  r  -  -T  '.uc 

MrO 

00  550  J=1,NYR3 
00  530  1=1, HOUR 

*'  =  V,  1 

IF  '  J f M , K ?  ,!.T.(-.l))  GO  TO  570 

X^Cr.O  CO  >  =  1 

0R(M.K)=9LANK 

A'IVrCI,K)s'NYR(I,K)y1« 

c  replace  flow  array  with  log  array 

TE^prJLCG (0(N,K) ♦00(1, K) ) *.4342345 
If  (IC0RL.E0.1)  3(H,K)=TEMP 
c  3UMt  SQUARES  AnO  CU3ES 

AV(I»N)3AV(I,K)*TE1‘P 
ro  CI,>0  =  SC  (l,n)*TE!'.P*TEOP 
o.nE  ia,X)=SKENCI,K)*TE«P*TErP*TEMP 
GO  SCO 

C  MIS  '  .  :W3  EQUATED  TO  T 

570  0(M,K).-T 
r,.?;.-,x)=E 
K^CP0=0 
580  CONTINUE 
59o  CONTINUE 

su*=o. 

00  t»?3  1  =  1  .NOUR 
•;“»=« LGG(I,K) 

IF  fTEMP.LT. 0.5)  GO  TO  820 
Tvr  sAVf :,X) 

AV(T,K)=TMP/TENP 

tr  C30(I, K) .LE.0.0.0R»TE1P.LT.2.5)  GO  TO  60C 
TVPA=30(I, K) 

s;)n,^)  =  (ona,K)-Av(I,K)*TNp)/{TEMP-l.7 
IF  (r.O(!,K)  ,Lt.O.)  GO  TO  600 
50(I,K)=SL'(I.K)**.5 


1207 

i2oe 

1209 

1210 
1211 
1212 
1213 
121« 

1215 

1216 

1217 

1218 

1219 

1220 
1221 
1222 

1223 

1224 

1225 

1226 

1227 

1228 

1229 

1230 

1231 

1232 

1233 

1234 

1235 

1236 

1237 

1238 

1239 

1240 

1241 
l?a? 

1243 

1244 

1245 

1246 

1247 

1248 

1249 

1250 

1251 

1252 

1253 

1254 

1255 
1 25b 

1257 

1258 

1259 

1260 
1261 
1262 

1263 

1264 

1265 

1266 

1267 

1268 

1269 

1270 

1271 

1272 

1273 

1274 

1275 
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SKEn(I,K)s.(TEMP*TEMp*sKEHa,K)-3,*TEMp*TMp*TNpA*2.*TMp*TKP*THP)/  1276 
TEMP-1. )*(TEHP«2.)  *30  CI,K)**3)  1277 

C"  TO  610  1278 

600  S0CI.K)=0.  1279 

SKSw(I,K)ao.  12S0 

610  SUM=oUM*5KEW(I»K)  128! 

620  CONTINUE  1282 

TEH2=N0UR  1283 

SU’^SUfVTEMP  128« 

N-N+l  1235 

IF(K.LT.NSTXX.ANO.NYOIP.EQ.O)  GO  TO  7 10  1286 

IF  (N.GT.l)GO  TO  630  1287 

C  PRINT  FREQUENCY  STATISTICS  1283 

WRITE (6, 1070) ISTA (K), (Av (I,K), Ial,NDUR)  1289 

WRITE (6, 1080) <SD(I,K),I=1,NDUR)  1290 

WRITE (6, 1090) (SKEW (I, K) ,1=1. NOUR)  1291 

WClTE(6rilOO)  (DQ(I,K),Tsl,NDtJ»)  1292 

WRITE (6, I l 10) C AN YU ( I, K). 1=1, NOUR)  1293 

IF  (ISWEVi.LE.Q. OR , INC A  0 .EE  .  0 )  GO  TO  710  1299 

630  IFCN.GE.16)  GO  TO  710  1295 

IF(SUM.GT.(AvG$K-.n.AN0.SUM.LT.(AVG3K*.l))  GO  TO  710  1296 

IN0C=1  1297 

M=0  1298 

DO  660  J=1,NYRS  1299 

00  650  1=1, HOUR  1300 

K“M*1  1301 

IF(u(N,K) .5E.T)  GO  TO  690  1302 

TcMP=Q(M,K)  1303 

0(M,K)=10.*«tEMP-0R(I,K)  1309 

GO  TO  650  1305 

690  Q(M,K)=-1.  1306 

650  CONTINUE  1307 

660  CONTINUE  1308 

IF  CS'JM-AVGSK)  670,710,690  1309 

670  DO  b50  1=1, NOUR  1310 

IF(NUOG(I,K).LE.O)  GO  TO  680  1311 

0Q(I,K)s0c(I,K)*1,5  1312 

680  CONTINUE  1313 

GO  TO  550  i j i q 

690  00  700  1=1, NOUR  1315 

IF  (nUCGCI,K).LE,0)  GO  to  700  131S 

OnCI,<)  =  0Qfi,K)-xiNCRCIrW)  1 3 1 7 

700  CONTINUE  1318 

GC  TO  560  U19 

710  CONTI N U C  1320 

IF(NYOIF.GT.O)  NSTXX=l  1321 

NSTAi  =  NSTA*Nr,TA  1322 

IF  (NOUR.EQ.l)  NSTAXPHSTA  _  _  ...  1323 

c  omit  correlations  if  only  i  sta  and  i  duration  i32« 

ITFNSSO  1325 

IF(ICORL.EO.l)  GO  TO  730  1326 

MsO  1327 

ANYR8=0.  1328 

DO  72C  J31 , NYR3  1329 

MrM+l  1330 

IF  CQCJfD.GE.T)  GO  TO  720  1331 

A  N  Y  .0  S  =  A  N  Y  (;  S  ♦  1  *  1332 

ISPI'-.D-sSLAnK  1333 

IRCKO (N) =1  1339 

720  CONTINUE  1335 

GO  TO  1760  1336 

C  OMIT  CORRELATIONS  IF  Nq  MISSING  FLOWS  1337 

730  IFfKRCRO.EQ.l)  GO  TO  1130  1338 

C  ******  COMPUTE  SUMS  OF  SQUARES  AND  CROSS  PRODUCTS  *******  1339 
790  DO  760  K=1,NSTA  j340 

00  150  1=1, NOUR  j3ai 

00  730  L=l,NSTAx  j3a2 

WA(I,K,L)=-«.  13«3 

SUMA ( I , K , L) =0,  i 344 

r.UM3(I.K,L)=0,  1345 
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S !3 A  f  I,  t; f  L)  =0. 

Cf  rI,K,L)=0. 

XPA3(I,K,L)=0. 

NCA9(I,K,L)=0 
750  CONTINUE 
760  CONTINUE 

00  OOC  ;<  =  1,MST* 

KX=K*1 

IF(KX.GT.NSTAX)  GO  TO  820 

M  =  0 

00  810  J=1,NYRS 
DO  800  I*1,NDUR 
M  =  n+1 

TE«P=3CM.K) 

IF (TEMR.GE.T)GO  TO  800 

I"  (ITHNS.EO.l)  TEMPsALCGdEMP*0Qd»>O)*.  *342945 
DO  790  L=X,NSTAX 

C  SURSCrIPTS  EXCEEDING  NSTA  RELATE  TO  AOJACENT  DURATION 

rra.LE.N5TM  SO  TO  770 
!.:<=L-nsta 

If  (I.EO.n  TMRrQ(H<-l  ,LX) 

IF£I.GT.i)TMP=GCM-t,LX) 
tr  (TMF.G'E.TIGO  TO  790 

:ra*!-:N3.r.Q.l)  7MP3ALCGdNP*0Qd» LX) 5  *.4342945 

GO  TO  7  80 

770  T'!R  =  0JN,L) 

•F(TNo,GE.T)GO  TO  790 

IFf.  JTRN5.E0.1)  THpsALOGITMp^Dod.L))*. 4342945 
C  COUNT  AND  USE  CNLY  RECOROEO  PAIRS 

730  N;4B(I,K,L)sNCAtl(I,K,L)*l 
3UMA (I ,K,L)SSbMA(IfK,L)7TEMP 

3UH3f  I,K,L)=9iJMB  CI,K,L)7TMP 
33  A  d,K/L)sSS«  (.tfK,U*TEMp*TEHP 
233  (I,K,L)=Sn*  (r,K,L)7TMP<*T«P 
xpaoc:,k,l'=xpa3{i,k,l)7TEmp*tmp 
IFa.GT.N3TA)  GO  TO  790 
NCARd,  L,K)=nCAB(I,k,L) 

3!.|MA(I.LrK)=Sl!NB{X,KfL) 

S  “bfT,L,  K)=SU*<Ad»K,L) 

:i'  ci.l.m=sob  (I,k.l) 

3/3  (la.X)=?CA  CI,K,U 
SCAB  ;  ►N)‘5xFAB  (*,K,L) 

790  C  r.TU'.E 

300  c;nt,u‘: 

e:c  cr.'vr *••'.£ 

c  »  *  *  *  *  «  COMPUTE  CORRELATION  COEFFICIENTS  ********** 

T  77PZ-J 

620  00  <’90  1*1 .  HOUR 

C  SEARCH  FOP  DURATION  WITH  LONGEST  RECORD 

I T  E  BPsslCO d » X) 

IFdTEMP.LE.ITMP)  GC  TO  830 

I rMP=iTENP 

IX=I 

830  ir (KX.GT.NSTAX)  GO  TO  870 
DO  360  L=KX,N3TAX 

C  ELTfiNATc  PAIRS  WITH  LESS  THAN  3  YR3  DATA 

TrfNCAR(I,K,L3 ,LE,2)  GC  TO  840 
Tc.PP  =  NCA6(I,KrL) 

SASSUMA(I,K,L) 

SR  =  S;JH3  ( I  r  K  .  L) 

)  :'p-(S0A(l,K,L)-3A**2/TEMP)*(S56df  K,L)-3B**2/TEHP) 

IF  £Tmp.LE.O.)  GO  TO  850 
THPB-1  . 

T"Pa  =  XPA9(I,K,U-3A*SB/TEHP 
IF(TRPA.LT.‘1.)TKP0s-TMPB 
T,|3A  =  Tf*PA»TMPA/TMP 

TMpA=l .-(l.-T, MPA)* (TEMP-1. )/(TEmP-2.) 

IF£TNPA.LT.0.)THPA=0. 

RA  d.A  ,L)=TMP£i*TMpA**.5 
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eao 

IN  CL.UT.NSTA)  Qj 0  TO  860 

1415 

ra(i,l,k)=sac:,k,l) 

1416 

GO  TO  860 

1417 

850 

RA(I,K,L)sO. 

1418 

860 

CONTINUE 

1419 

c 

ELIMINATE  NEGATIVE 

GROSS 

CORRELATIONS 

1420 

870 

DO  880  L=1,NSTAX 

1421 

1  •'  t  >  3 

if  (tekp.lt. o.o. and. temp. ge 

.(-1.0))  RA(I,K,L)=0. 

1423 

800 

CONTINUE 

1424 

BA(I,K,K)sl. 

1425 

890 

CONTINUE 

1426 

900 

CONTINUE 

1427 

IF  ( I TRNS • NE  .0)  GO  TO  1270 

1428 

c 

•  * 

*  •  *  «  ADJUSTMENT  OF  FREQUENCY 

STATISTICS  TO  LONG  TERM 

1429 

DO  980  II=1,NDUR 

1430 

I=IX+II-1 

1431 

IFCI.GT.NDUP) I=N0UR-II*1 

1  ®32 

00  910  K=1,N3TA 

1433 

ISTN(K)SX 

1434 

ISTy (K)*NL0G(I,K) 

1435 

910 

CONTINUE 

1436 

c 

ARRAY  STATIONS  -  LONGEST 

RECORD  FIRST.ETC 

1437 

I7*P=NSTA-i 

1438 

IF  (ITMP.LE.0)  GO  TO  985 

14361 

DO  930  KX=1,ITMP 

1439 

ITP=<X+1 

1«40 

DO  920  K=ITP,NSTA 

l«4l 

IF(ISTY(KX).GT.ISTV(K))  GO 

TO  920 

1442 

ITEMP=ISTU(KX) 

1443 

ISTNCKX)=ISTn(K) 

1444 

I3TN(K)=ITEMP 

1445 

ITENP=rSTY(KX) 

1446 

ISTY(KX)=ISTY(K) 

1447 

ISTY(K)=ITEMP 

1446 

920 

CONTINUE 

1449 

930 

CONTINUE 

1450 

on  970  KX=1,N3TA 

1451 

K=ISTN(KX) 

1452 

TMP3=NLCG(I,K) 

1453 

IN0C=0 

.  ... _ 

1454 

DO  960  LX=1,KX 

1455 

IF  (LX.EQ.KX)  GO  TO  940 

1456 

ITP=I 

1457 

L=ISTN(LX) 

1458 

Tmp=MLCG(I,L) 

1459 

TMPP=NC*B(I,k#L) 

1460 

GO  TO  950 

1461 

940 

IF(NDUR.EQ.l)  GO  TO  960 

1462 

ITP=I-1 

1463 

IF(ITP.LE.O)ITP=I+l 

1464 

L=K+NSTA 

1465 

TMP=NLOG(ITP,K) 

1466 

TMPP=NCA8(I,K#L) 

1467 

950 

TP=RA ( I , K , L) 

1468 

IFCTP.LT. (-1.))  GO  TO  960 

1469 

TMPA  =  TMPP/n.-(THP-TMPP)«TP»*2/THP) 

1"70 

IFCTMPA.LT. TMPB)  GO  TO  960 

1471 

INDC=1 

1472 

ANYR(I,K)=TMPA 

1473 

TMPBsTMPA 

1474 

ITMprU 

1475 

ITEMPSITP 

1076 

960 

CONTINUE 

1477 

IFCINOC.LE.O)  GO  TO  970 

1478 

LsITMP 

1479 

ITPsITEMP 

1480 

LX=L 

1461 

IF(LX.GT.NCTA)  LX=LX-NSTA 

1462 

TP  =  (.;A  (I,  K,L) 

1063 

TE9Pa,NCAP(l,K,L) 

1484 
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S'.-'U'V*  ll,K,L) 

S«  =  SL‘MP  (I,K,L) 

S0A=fsr;A(I,K.t)-3A*«3/TEMP)/(TEMP-l.) 

IF(SDA.LT.O.)  SDArO. 

S0A=30A«*.5 

S00=(SQ8(  I,K,U-SB*«2/TEMP)/(TEHP-1.) 

IF  (SD8.LT.  .0005)  GO  TO  <370 

S0?=SDR**.5 

THPP-S0A/3DB 

AV(I»K)=SA/TFMP+(AV(lTP»LX)»3B/TEHP)*TP*TMPP 
SO(I,K)=SDA*(SO(ITP,LX)-SOB)*TP**2«TMPP 
970  CONTINUE 
980  CONTINUE 

985  IF  (ISKEW.GT.O)  GO  TO  1020 

lFCNSMTH.LE.(-l))  GO  TO  1050 

SMOOTH  3KEW  COEFFICIENT 
DO  1080  Ksl , N3T A 
SA  =  0. 

3C  =  0. 

SAA=0. 

SAC  =  0. 

I TMp=NDUR 
DO  1000  1=1, NDUR 
IF(NLCG(I,tO.LT.3)  GO  TC  990 
IF(S*Ew(I,K).GT.l.)  SKEw(I,K)=l. 

IF (S«EhCI,K) .LT.(-1.)J  3KEw(I»K)a-l. 

IF  (  NDUR . LT . 3)  GO  TO  100Q 

tp=av(i,k)-auoG(p(i)) 

TEMP=SKEH(I,K) 

SA  =  SA-fTP 
SC=SC+TEMP 

sa*=saa+tp«tp 
sac=sac+tp*temp 
GO  Tn  1000 
990  ITfP=ITMP-l 
1000  CONTINUE 

IF(IT MP.LT, 33  GO  70  1050 

TP- I TMP 

SAA=oAA-SA*SA/TP 
SAC=3AC-5A*SC/TP 
BC=3AC/S«A 
IFOC.GT.l.)  BC  =  1. 

IFf8C.LT.  (-1.3)  8C«-1. 

CC*(SC“BC*SA)/TP 
DO  1010  1=1, NDUR 
TEMP=AV(I,K)-ALOGCP(X)) 

SKEN(I,K)=CC+eC*TEMP 
1010  CONTINUE 
1040  CONTINUE 

GO  TO  1050 

1020  DO  1030  I=1,NDUR 
DO  1030  K=1,nSTA 
SKEW(I,K)=SKWU) 

1030  CONTINUE 
1050  WP7TSC6, 1060) 

1060  FORMAT  C/63H  FREQUENCY  STATISTICS  after  ADJUSTMENT  WITH  A  LONG  term 

I  STATION  ) 

HR  HE (6, 990) (A A (I) ,A8(I),I=l,N0uR) 

DO  1120  K=1,n3TA 

W9lTE(6,ln70)ISTA(K), ( A v ( I , K ) , 1= 1 , NDUR) 

1070  FORMAT (/ 10, BH  MEAN  1CF11.3) 

WRITE (6,  1080)  (SO  (  I, K 3,1  =  lr NDUR ) 

1080  F0RMATf9x,7HST0  DEV  10F11.3) 

WRITE  (6, 1090)  (SKEW  f I , K) ,1  =  1, NDUR) 

1090  F0R'1AT(12x,4HSKEW  10F11.3) 

KPTTE(6, 1 100) (00(I,K) , 1=1 ,NDU«) 

1100  FORMAT (10x,6HlNCRMT  F10.2,9F11.2) 

WRITE (6, 2000) (ANYR(I,K),I=1,NDUR) 

1110  Format  (Ux,5Hyear3  lOFll.O) 

DO  1120  1=1, NDUR 
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AufR(I.K)=NLCG(I,K)  1559 

t 1 20  CONTINUE  1555 

1556 

c  ******  TRANSFORM  TO  STANDARDIZED  variates  ***********  1557 
1130  00  1160  K=1,NSTA  1558 

MrQ  1559 

DO  1170  J=1,NYRS  1 5b0 

CO  1160  I=1,NDUR  1561 

Msm+1  1562 

IFCO(M,K).GE.T)GQ  TO  1160  1563 

IF(sn(I,K).i_E.0.)G0  TO  1150  1569 

Q(M,K)a(Q(M,K)-AV(I,K))/SO(I,K)  1565 

C  PEARSON  TYPE  III  TRANSFORM  1566 

TMPPs  SKE  W ( I »  K )  1567 

IF (TMPP.EQ.O.)  GO  TO  1160  1566 

TEMP=.5*TMPP*R(M,K1+1.  1569 

TKP=1.  1570 

IFfTEMP.GE.O.lGC  TO  119C  1571 

TEMP=-TEMP  1572 

TNPs-T"P  1573 

1 1  AO  Q(“,K)=6.*(TMP*TEmp*«(1./3.)-1.)/TMPP  +  TMPP/6.  1579 

GO  TO  1160  1575 

1150  0(MrK)=0.  1576 

1160  CONTINUE  1577 

1170  CONTINUE  1578 

1160  CONTINUE  1579 

ITRNS=-1  1560 

GO  TO  790  1561 

C  ******  ESTIMATE  MISSING  CORRELATION  COEFFICIENTS  •  **•***»  1562 
1190  IF(NSTA.LE.l)  GO  TO  137C  1583 

00  1260  I  =  l,NDlJR  1569 

IX=:-1  .  1585 

IF(I,EQ.1)IX=I*1  1566 

DO  1 P50  K=1,N3TA  1587 

K  X*K ♦  1  1 56B 

IF  (KX.GT.NSTAY)  GO  TO  1250  1569 

00  1290  L=KX,MSTAX  1590 

C  L  anO  K  CORRELATION  POSSIBLY  MISSING  1591 

Ie(RMI,K,L).GE.<-l.))GC  TO  12«0  1592 

RMAX=1.  1593 

RMIN=-1,  1599 

C  LX  SEARCHES  ALL  OIRECTLY  RELATED  CORRELATIONS  1595 

DO  1230  LX=1,NSTAX  1596 

IF(LX.EQ.K)GO  TO  1230  1597 

IFfLX.EO.UGO  TO  1230  1598 

TtMP=RA(I,!<rLX)  1599 

IF (L.LE.NSTA)GO  TO  1200  1600 

IF(LX.LE.NSTA)SC  TO  1210  1601 

c  both  l  ano  lx  represent  adjacent  durations  1602 

ITMP  =  L-NS7A  If  03 

ITEMP=LX-NSTA  1609 

TMPSRACIX, ITMP. I  TEMP)  1605 

GO  TO  1220  1606 

c  l  represents  current  duration  1607 

1200  TMP=RA(I,l,LX)  1606 

GO  TO  1220  1609 

C  LX  ANO  NOT  L  REPRESENTS  CURRENT  DURATION  1610 

1210  TKP=?A{! rLX,L)  1611 

1220  IF (tmp*tEmp.lT. (-2.) )G0  TO  1230  1612 

TMPA=( (1 .-TEMP»TEMP)«(1 ,-TMP«TMp))»*,5  1613 

T';pp,itmp«temp*tmpa  16io 

:FCTwP9.LT.HMAXJRMAX3TMpB  1615 

Tr:Pt)  =  Tf'Pe-TMPA-TMPA  1616 

IF(TMPS.GT.PMIN)9'1IN  =  TMP8  1617 

1230  CONTINUE  1616 

C  AvFRAGE  SMALLEST  MAX  AND  LARGEST  MIN  CONSISTENT  VALUE  1619 

RACI»K,L)  =  CRSAX*Rr1IN).,5  1620 

IF  (RA(I,K,L)  .LT.O.O)  R*(I,K,L)»0.  1621 

IF(L.LE.N3TA)KA(I,L,K)SPACI,K,L)  1622 
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1240  r-.Mt'-'J- 

1250  CONTINUE 
1260  CONTINUE 

GO  TO  1 3 T 0 

C  »«***«  PRINT  CORRELATION  MATRIX  **************** 

1270  DO  1360  I  s  1 , NDUR 

IE(ITRNS.LT.l)  WRITE (6, 1280) A A (I) ,AB(I) 

1280  FOrma T ( //46H  CORRELATION  COEFFICIENTS  OF  RECORDED  OA+A  FOR  2A4,9H 

1  CIJR  NT  I  ON  ) 

lr  (  ITKNS.GT.01  WRITE (G, 1290) AA  (I) , A8(I) 

1290  f.irnat ;//6oh  correlation  coefficients  of  recorded  and  reconstitute 

ID  DATA  FOR  2A4.9H  DURATION) 

WRITE (6#  1300)  CISTA  CK) ,K=1 , N3TA) 

1300  FORMAT  (/3X,3HSTA  2016) 

WRITE (6, 1310) 

1310  FORMAT  (20x,18HwITH  5  *  ME  DURATION) 

00  1320  K=1,NSTA 

WRITE (6,1 330) ISTA(K), (RA (I,K,L) ,L=l,NgTA) 

1320  Continue 

1330  Format Cix, I5.20F6.3) 

IF  (NOUR.EQ.l)  go  TO  1360 
WRITE  (6, 1340) 

1340  Format (20X,39HWITH  ADJACENT  OURATION  AT  ABOVE  STATION) 

ITFsNST  A* l 
DO  1350  K=i,nSTA 

WRITE  (6, 1  330)  1ST* (K) , (RA (I,K,L) ,L  =  ITP,NSTAX) 

1350  CONTINUE 
1360  CONTINUE 

IF  (HRCbQ.E3.1J  GO  TO  1760 
IrCTRNS)  1  190,1190,2020 

C  *«,»«*  RECONSTITUTE  MISSING  OATA  *••*••**•••*•••• 
1370  M=0 

Nvar=NSTA+i 

DO  1750  J=t,NYRS 

IF  (IRCRO(J).EQ.l)  GO  TO  1380 

m=m+nouR 

GO  TO  1750 

1380  DC  1740  1  =  1  ,N'OUR 
mm-m 
MX  =  M 
M  =  M*1 

IF  CT.E3.1) MMaM*l 
On  r/30  K=1,NSTA 
KX=NSTA*K 

IF  (i)(M,K)  .LT.T.OR.NLOG(I,K)  .LT.3)  GO  TO  1730 

Nir;op  =  0 

I PRE  VsO 

C  form  CORRELATION  MATRIX  FOR  EACH  MISSING  FLOW 

OD  1450  L=1,NSTA 
LA  =  MI  NOP 

IF  (L.cO.K)GO  TO  1420 
IFfQ(M.L) .GE.T)GO  TO  1450 
NINDP=NINDP*1 
X(Nf  :JP)=rj(M,L) 

DO  mo  LX  =  L ,  NST  A 
IF (LX.EQ.K)GO  TO  1390 
IF(0(“,LX),GE.7)GO  TC  1410 
LA=l»*1 

R(NINDP,LA)=rA(I,L,LX) 

GO  TO  1400 

1390  IF  (NDUR.E0.1J  GO  TO  1410 

IF(R(MM,LX).GE.T)GO  TO  1410 
la=la*i 

R(NI‘0°,LA)  r  RA(I,L,KX) 

1400  R(LA,MNOP)  e  R ( NI NDP » L A  ) 

1410  CONTINUE 

R(NINDP,NVAR)SRA(I,L,K) 

GO  TO  1450 

1420  jc  (MDUR.EQ.I)  GO  TO  1450 
IFCDfMM.M) .GE.TJGO  TO  1450 
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NINDP=MIN0P*1 
IP»tV=NINDP 
X  (NINDP)=U(MM,L) 

00  1440  LX  =  Lf N3TA 
IP  (LX.EQ.K)GO  TO  1430 
IF(a(M,LX)  .GE.TJGO  TO  1440 
la=la+i 

R(NINOP.CA)  =  RA(I,LX,KX) 

R(LA,NINDp)=R(NINDP,LA) 

GO  TO  1440 
1430  LA=IA+1 

R(NINOP,LA)sl . 

1440  CONTINUE 

R(NInDP,NVAR)*RA(I,L»KX) 

1450  CONTINUE 

C  CASE  NUMBER  1  RESULTS  WHEN  NO  FLOWS  ARE  FOUND  FOR  CORREL* 

ICSE=1 

IF(nINDP.LE.O)  GO  TO  1510 
ITMP=NINOP-H 
DO  1460  IXsl.NINDP 
1460  R(IX, ITMP)=R(IX,NVAR) 
c  ========== 

1470  CALL  CRDUT(R) 

C  ========== 

ITE“P=NINDP+1  _  _ 

TEMP=1. 

IN0C=0 

00  1 4 p 0  L  =  1,NIN0P 
T  H  p  =  A  B  S  (R  CLr  ITEMP)  ) 

IF(TMP.GT.TFMP)  GO  TO  1480 
:F(L.EQ,lPREV.AN0.TMP.GE..9)  GO  TO  14S0 
TEMPsTMP 
ITP  =  U 

1480  Ir(R(L,ITEMP).LT.0..AND,8(L).6T,(«l»5).AND.B(L).LT».5)  GO  TO  1490 
IF  'R'L.ITEMP).GT.0,.AN0.9(U.GT.(-.5).AND#B(LJ.LT.1.5>  GO  TO  1490 

Innc-l 

1090  CONTINUE 

IFfINDC.GT.OJ  GO  TO  150C 

IF(DTR«C.LE.1..AND.0TRMC.GE.0.)  GO  TO  1590 
C  IF  matrix  INCONSISTENT,  0HIT  VARIABLE  WITH  LEAST  CORRELAT 

1500  IT“P=NIN0O-l 

IFflTMP.GT.O)  GO  TO  1530 

C  CASE  NUMBER  2  RESULTS  WHEN  ALL  CORRELATIONS  ARE  ZERO 

ICSE=2 

c  possible  branch  from  870*2 

1510  IyR=IYRA+J 

WRITE  Cf>, 1520)  ISTA  (K)  ,  I#  I YR,  ICSE 

1520  FORMAT (/25H  ZERO  CORRELATION  FOR  STA  ,I6,10H  DURATION  ,12, 6H  YEA 
1R  •  1 5 p  6 H  CASE  , 12/ J 
PC1)=0. 

X ( ! ) -0  . 

0TRMC=0. 

GO  TO  1590 

1530  IFUTP.3T.ITMP)  GO  TO  1580 
DO  : 55C  L=!TP,ITMP 

oi  l5uo  la=i,;tem? 

1540  P(L,LO=R(L  +  l,LA) 

1550  r(L)=X(LM) 

1560  OH  1560  L=1,ITMP 

00  1570  LA=ITP,NINDP 
1570  RCL,LA)=R(l,LA+1) 

1580  CONTINUE 

NINO.=>=ITMP 
GO  TC  1470 

C  ADO  RANOOM  COMPONENT  TO  PRESERVE  VARIANCE 

1590  Tmo=RNGEN(Ixx) 

■’•EMrrRsGSN(IxX) 

TEMP=(-2.*aL0G(TEMP))**,5*3IN(6,2832*THP) 

C  COMPUTE  FLOW 
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EXHIBIT  6 


T£MP=TEMP«C1  .-PTRMC)**.5 
DD  1 S00  L=t,NlNDP 
TEHP=TEHFt3{L)*X(L) 

1600  CONTINUE 

0(M,k)=T£MP 

ANYR(I,K)sANYR(I,K)+DTRMC 

TP=Q(M,X) 

C  ADD  NEW  VALUE  TO  SUMS  OF  SQUARES  AND  CROSS  PRODUCTS 

00  1670  L=1,NSTAX 

C  SUBSCRIPTS  exceeding  NSTA  relate  to  preceding  month 

1610  IF(L.LE.N3TA)  GO  TO  162C 
LX=L-NSTA 

IF  (I.EQ.l)  THP=Q(M*1,LX) 

IF(I.GT.l)  TMP=Q(M-1,LX) 

GO  TO  1630 
1620  TMP=Q(M,L) 

1630  ir(TMP.GE.T)  GO  TO  1670 
C  COUNT  AND  USE  CNLY  RECORDED  PAIRS 

NCA8(If  k,U=nCAB(I,K,L)  +  1 
SJMACI,KfU=SUYA(I,K,L5+TP 
S'JM3(I,K,L)=SUM3(I,K,L)«TMP 
SQA  (I,K,l5=SG*  CI,K,U+TP*TP 
S'.’B  CI,K,L)=S08  ci,k,l)  ♦thp*tkp 
XPAB(I,K»L)=XPA3  (  I «  K  » L  3 ♦  TP*  T  HP 
if  CL. GT. NSTA)  C-0  TO  1640 
NCAB(I,L»K)snCAB(I,K#L) 

SIJNA(I,L,K)=SUMB(I,K,L) 

S’JMB  (I,L»K3=SUMA(I,K,L) 

S3A  (I,L,K)=SUB  (I,K,L3 
SOB  (1,L,K)=SUA  (I,K,L) 

XPAB(I,L,K)=xPAB(I,K,L) 

C  RECOMPUTE  CCRRELATION  COEFFICIENTS  TO  INCLUDE  NEH  DATA 

C  ELIMINATE  PAIRS  WITH  LESS  THAN  3  YRS  DATA 

i 6 AO  IF(nCAB(I,K,L).LE.2)  GO  TO  1670 
TEYr  =  NCAOf.Ir<,L) 

TMPs(SQA(I,K,l  ’  ..  .  ,YtKfL)*SUMA(I,K»L)/TEMP)»'<*QB(I,K»L)-SUMB 

1 ( I »  K , L  3  *  SUMR 1  '  ’  '  £Mp ) 

C  ELIMINATE,  pairs  with  ZERO  VARIANCE  »  ,<luuCT 

IF  (TMS.LE.O.)  GO  TO  lfc 50 
TUP3=1 . 

tmpa  =  xpaB(I,k,U-SUMA(  t.K,L)«SUM  I,KfL)/TEHP 
c  retain  a l g e r r  *ic  sign 

IF  (TMPA.LT. n.)TMPEs-TF  J8 
T(1PA  =  TMPA*TU0A/TMP 
RA(I,K,L)=TMpB«TMPA**,j 
IF(RA(I,K,L).GE.0.)  GP  to  1660 
1650  RAfI,K,L3=0. 

1660  IF (L.GT.NSTA)  GO  TO  1670 
RAfI.L,R)sRA(I,K,l ) 

1670  CONTINUE 

IF(NDUR.EQ.1)G0  TO  1730 
00  1720  L=1,NSTA 
TTP  =  0 

I x- 1 ♦  1 

IFCIX.GT.NDIJR)  GO  TO  16e0 
TMr=Q(M*l,L) 

GO  TO  1700 

1680  IC(I.GT.2)  GO  TO  1730 
1690  TNP=0(M-1 ,L) 

IX=I-1 

ITP=1 

1700  IF  (T^P.GE.T)  GO  TO  1720 

NCAH(Ix,L,KX3cNCAB(IX,L,KX)*l 
SllMA  (IX,L,KX3=SIJMA(IX,L,KX)*TMP 
SUMS (IX, L,KX)=SUMBCIX,L,KX)YTP 
SQA  (IX, L, KX)=SQA  (IX,L,KX)*TMP**2 
S'JB  (IX,L,KX)=SQB  (IX,L,KX)*TP**2 
XPAt3(Ix,L,KX)aXPAB(IX,L,XX3+TMP«TP 
IF(NCAB(IX,L,KX).LE.2)  GO  TO  1720 
TEMP=NCA9(IX,L,XX) 


I7fc2 

1763 

1764 

1765 

1766 

1767 
1760 

1769 

1770 

1771 

1772 

1773 

1774 

1775 

1776 

1777 

1778 

1779 

1780 

1781 
t’32 
1783 
1764 

1785 

1786 

1787 
1786 

1789 

1790 

1791 

1792 

1793 

1794 

1795 

1796 

1797 

1798 

1799 

1800 
1801 
1802 

1803 

1804 

1805 

1806 

1807 

1808 
1  809 
1810 
1611 
1012 
1613 

1314 

1815 

1816 

1817 

1818 

1819 

1820 
1821 
1822 
1823 
1624 

1825 

1826 
1827 
1626 

1629 

1630 

1631 
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EXHIBIT  fi 


TMP=CsQl(IX,UfKx)-SlJHACIX,L,KX)**2/TEMP)*(3QB(IX,L»KX)-  1«32 

lSUM3(IX,L,KX)**2/TEi1P)  1633 

IF(TMP.LE.O.)  60  TO  1710  183'' 

TMPB=1.  1635 

TrlPAsxPAB(IX.L,KX)-SUMA(IX,L,KX)*3UMB(IX,i.»KX)/TEMP  1336 

IFCTMPA.LT.O.)  TMPB3-TMP8  1637 

TMPA=TMPA**3/TMP  1638 

RA(IX,L»KX)=TMPB«TMPA**.S  1639 

IF(RA(IX,L,KX).GE.O.)  60  TO  1720  1640 

1710  RA(IX,L,KX)=0.  1841 

IFCI.E0.2.AN0.ITP.LT.1)  GO  TO  1690  1842 

1720  CONTINUE  1843 

1730  CONTINUE  1844 

1740  CONTINUE  1845 

1750  CONTINUE  1846 

1760  WSItE(6,50)  1847 

WSITE(6, 1770)  1648 

1770  F0RMATI33H  RECORDED  AND  RECONSTITUTED  DATA  .  1849 

DO  1930  K=1,N3TA  1850 

IFCX.GE.NSTXX)  WRITE(6,  1780HAA(I),AB(I),I*1,NDUR)  1851 

1780  FQRMAT(/2X, 10H  STA  YEAR  4X , A3 , A4 , 9  ( 3x , 2 A4 ))  1852 

M  =  0  1853 

C  CONVERT  STANDARD  DEVIATES  TO  FLOWS  1854 

ANYR3=NYRS  1855 

00  1890  J=1,NYRS  1856 

IF  (IRCRO(J)  .Efl.l)  GO  TO  1790  1857 

M  =  M«.NDUR  1858 

ANYR3=ANYRS-1.  1859 

GO  TO  189o  I860 

1790  DO  1870  I=1,N0UR  1861 

M=M»l  1862 

XCn=OR(H,K)  1863 

XOm=a(MfK)  1664 

IF (ICURL.EQ.OlGO  TO  1870  1865 

IF  (NLOGCI.K) .LT.3)  GO  TO  I860  1666 

TEMPnO(NfK)  1867 

TMP=SK£W(I,K)  1668 

C  USE  ADOPTED  SKEW  FOR  RECONSTITUTING  1869 

IF (I  SKEW. GT. 9)  TMP  =  SKW(n  1870 

IF(TMP.E3.0.1  GO  TC  182C  1871 

TF.MP=(  (T*p*  (tEmP-TmP/6,)/6.'»1.)**3-1.)*2./TMP  187  2 

IF (ORCP»K) .NF.E)  GO  TO  1820  1873 

T«pp=(-2,)/T*P  1874 

IF(TMP)  1800,1820,1010  1875 

1800  IF(TEMP.GT.TmPP)  TEKPsTPPP  1876 

GO  TJ  1820  1877 

1810  TF (TEMP.LT.TNPP)  TEMP=TPPP  1878 

1820  TMP=TEHP*SD(I»K)4Av(IfK)  1879 

TEPP=10.*»TPP-D  1(I#K)  iaeo 

IF (TE'P.LT.O.)  TEMP=0,  1881 

IFITENP.LT. QMIN(I,K))  QKIN(I,K)STEMP  1882 

Q(H,K)*TtMP  1863 

IF(I.EC.l)  GO  TO  1850  1884 

TMP  =  <11N-I,K)*P(I>/P(I-1)  1885 

IF(0(M,K).LT.THP)  GO  TO  1850  1886 

IF(QR(*,K).EC/.E)  GO  TO  1840  1867 

ITP=I-1  1888 

00  1 5 30  L  =  1,ITP  1869 

TNP=GCN-L,K)*P(I)/P(I-L)  1890 

IF(TNP.LT.O(M, K).AND.QR(M-L,K).EG.E)0(N-LfK)BO(H,K)*P(I-L)/P(I)  1891 
IF(ML0G(I-L,K).GT.2)  XQ(I”L)sQ(M*>L#K)  1892 

1830  CONTINUE  1893 

r-0  T'J  1850  1894 

1840  GC*,K)*THP  1895 

1850  XO(I)aG (M,K)  1896 

GO  TO  1870  1897 

i860  Xim=-1.  1898 

1870  CONTINUE  1899 

IF(K.LT.NSTXX)  GO  TO  1890  19C0 

IVR=IYRA*J  1901 
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EXHIRIT  6 


WRITE (6,  1*60)  I3TA(K)  ,IyR,  CXQCI)  f  X(I)  #  1=1, NDUR) 

1902 

1880  FORMAT (2 16, FI  1.0, Al.F 10.0, A  1,8 (F10. 0, Al)) 

1903 

Iff IPCHQ.GT.O)  WRITE (7, 60)  I  ST A (K) , I YR , ( XQ ( I) , 1  =  1 , NDUR) 

1904 

1890  CONTINUE 

1905 

IF(X.LT.NSTXX)  GO  TO  19e0 

1906 

IF (ICORL.EQ.O.OR.KRCRD.GE.l)  GO  TO  1980 

1907 

IN0C=0 

1906 

1900  DO  1 R 1 0  1  =  1, HOUR 

1909 

IF  (DM  I N  (I,K)+D3(I,K) .GT ..0001)  GO  TO  1910 

19 1 C 

INDC=1 

1 9 1 1 

1910  CONTINUE 

1912 

IF(INOC.LT.l)  GO  TO  1930 

1913 

DO  1920  1=1, NDUR 

1914 

1 

D3(I,K)=00(I,K)+XINCR(I,K) 

1915 

1920  CONTINUE 

1916 

GO  TO  1900 

1917 

*  »  1 9 1 8 

1930  DO  1970  I=1,NDUR 

1919 

IF  (NL0G(I,K).LT.3)  GO  TO  1960 

1920 

TMP=0. 

1921 

TEMP=0. 

1922 

THPA=0. 

1923 

M=I 

1924 

DO  1950  J=1,NYRS 

1925 

IF  (IRCRD(J) .EO.O)  GO  TO  1900 

1926 

TP=ALOG(Q(M,K)fDQ(I,K)) 

1927 

TmP=TMP*TP 

1923 

TEMP=TEMP*TP*TP 

1929 

tmpa=tmpa>tp*tp«tp 

1930 

1940  M  =  M  ♦  NDUR 

1931 

1950  CONTINUE 

1932 

AV(I,K)=TMP*. 4342945/ ANYRS 

1933 

S0(  !,*)  =  (  (TF.MP-TMP<*TMP/ANYRS)/(ANYRS-1.))**.  5 

1934 

SKEwC,K)=(ANYRS»ANYRS*TNPA-3.*ANYRS*TKP*TEMP+2.*TFP*»3)/ 

1935 

1  (ANYRS»(ANYR3-1.)*(ANYR3-2.)*S0(I,K)»*3) 

1936 

1937 

SD(I,K)=SD(I.K)*. 4342945 

1938 

C-0  TO  1970 

1939 

I960  ANYR(I,K)=0. 

1940 

1970  CONTINUE 

1941 

1900  CONTINUE 

1942 

IFCCCPL.EO.O.OR.KRCRD.GE.l)  GO  TO  2020 

1943 

wrITE(6,50) 

1944 

WRITE(6, 1990) 

1943 

1990  F Oft NAT (//56H  FREQUENCY  STATISTICS  OF  RECORDED  AND  RECONSTITUTED 

OA  1946 

1  T  A  ) 

1947 

WRITE  (6, 490) (AA(I),AB(I),I  =  1,N0UR) 

1940 

00  2010  KsNSTXX.NSTA 

1949 

WRITE (6, 107  0) I  ST A(K), (Ay(I,X),I  =  l, NDUR) 

1950 

WRITE  (6, 1080) ( SO ( I, K), 1=1, NDUR) 

1951 

WRITE  (6, 1090)  (SKEW(I,K),I  =  1,N0UR) 

1952 

WRITE  (6,2000)  (ANYR(I,K), 1=1, NDUR) 

1953 

2000  FORMAT  ( 7 x , 9HE0U I V  YR8  10FI1.1) 

1954 

2010  CONTINUE 

1955 

C  RECOMPUTE  CORRELATION  MATRIX 

1956 

ITRNS=1 

1957 

• 

GO  TO  730 

1958 

*  *  1959 

2020  IT"T'  =  At- YK3»  .1 

I960 

C  COMPUTE  MEDIAN  PLOTTING  POSITIONS 

1 9b  1 

• 

TEMr=l ./ANYR3 

1962 

PLTT(1)=(1.-.5»*TEMP)w100. 

1963 

TE^P=(1C0,-PLTT(1)-PLTT(1))/(ANyR3-1.) 

1964 

00  2030  J=2 , ITMP 

1965 

FLTT(J)=PlTT(J-1)»TEMP 

1 9b6 

2030  CONTINUE 

196  7 

WRITE (6,2040) 

1968 

2040  FORMAT (//17H  FREQUENCY  ARRAYS) 

1969 

DO  2130  X=NSTXX,N3TA 

1970 

DO  2080  1=1, NDUR 

1971 
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EXHIBIT  6 

M=t  1972 

Gm(I)=SCM,K)  1973 

IF  (QM(I) .GE.T)  QM(I)»-T  1974 

X(I)=Qfi(K,K)  1975 

JA=1  1976 

00  2070  J=2,NYRS  1977 

M  =  m  +  N0UR  .  »•  *  1978 

IF  (IHCRO(J).  .  J.  -  .1  TO  2070  *  ’  1979 

JA=JA+l  I960 

TEMF=0(M,K)  1961 

JX=JA*NDUR+I  Jf . 1  1982 

DO  2050  L=2,JA  19„3 

LX=JX-L*NDUR  1984 

ITP=LX*MOUR  1965 

IF  COM (LX) .GF. TEMP) GO  10  2060  1986 

O'l  Cl  TP) SQM (LX)  19«7 

X(!TP)=X(LX)  !9ae 

£050  CONTINUE  1Q69 

QN(I)=TEKP  l990 

X(I)=QR(H,K)  1991 

GO  TO  2070  1992 

2060  QM(ITP)=T£MP  1993 

X(ITP)=QR(M,K)  1994 

2070  CONTINUE  1995 

2000  CONTINUE  1996 

WR1T£(6,2410)ISTA(K)  1997 

2090  F0RMATC/10H  NO  PLOT  3X , A3 , A4 , 9 (3x , 2 A4) )  1998 

WRITE  (6,2090) C A A ( I ) , AD C I) , I* 1 , NOUR )  1999 

H=n  2000 

DO  2120  J=t,ITMP  2001 

DO  2100  I=1,NDUR  2002 

KSM+1  2Co3 

XCn=X(M)  2  0  0  4 

X0CI)S3M(M)  2005 

IF(NLCG(I,K).LT.3)  Xfl(I)«-l,  2006 

2100  CONTINUE  2007 

WRlTE(6,2H0)JtPLTT(J),  (XQ(I},XCI)f  JsirNOUR)  2008 

2U0  F0RMAT(1X,I3,F6.2,F11.0,AI,9CF10.CA1))  2009 

2120  CONTINUE  20io 

2130  CONTINUE  20-1 

GO  TO  219Q  2012 

C  ft**********************************  2013 
C  RE  AO  STATISTICS#  IF  SUPPLIED  2014 

2140  WRI TE  (6 »  2  ISO)  2015 

2150  F0RMATC/27H  INPUT  FREQUENCY  STATISTICS  )  2016 

WRITt(6,490)  (AA(I),AB(I), 1=1, NOUR)  20l7 

OD  2180  K=l,N3TA  2018 

DD  2170  1=1, NOUR  2019 

c  «*  CARD  I  **  2020 

READ  (5,2160)  XSTA (k) , A V ( I , K) , SO ( I , k) , SKEW ( I , K) » OQ ( I , K ) , A N YR ( I , K)  2021 
2160  F0R“ATClX,I7,8xf5F8.0)  2022 

NLOGCI,K)=ANYR(I,k)  2023 

2170  CONTINUE  2024 

WRI TE( 6, 1070) I  STACK), C A V C I , K ) , 1  =  1 , HOUR)  2025 

WRITE (6,1080) (30(1, K), 1=1, NDUR)  2026 

WRITE (6, 1090) (SKEw(I,K) ,1=1, NOUR)  2027 

WRITE  (6, 1100)  (t)Q(I,K)  ,1  =  1,  NOUR)  2028 

WRITE (6,2000)  (ANYR(I,K),I=1,N0UR)  2029 

2180  CONTINUE  2030 

2190  Dn  2250  K=n3TXX,N3TA  2031 

C  ******  SMOOTH  STATISTICS  ******•***•»##**#*»*  2032 
IF  (NSPTH.lE.C-1))  GO  TC  2230  2033 

IFCNDUR.LT. 3)G0  TO  2230  203« 

C  SUMS#  SQUARES  ANO  CROSS  PRODUCTS  2035 

St*0.  203b 

S0=0.  2037 

SC=0.  2038 

S  A  A  =  0  .  2039 

SA8=0.  2040 


SAC=0. 

ITMPsNOUR 


2041 

2042 
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EXHIRIT  6 


CO  2210  IaIfNDL'R 

2083 

IF  (NLCGCI.K) .LT.3)  GO  70  2200 

2088 

TP=AVCI,K)-ALOGCP(I)) 

2085 

TnP=SO(I,K) 

2086 

IF  (SKcXI.K)  .GT.l.)  SKEt,(I,K)  =  l. 

2087 

I- (?KEwCI,K).LT.(-l.))  SKEw(I,K)»(-l.) 

2048 

TlM3=SKEW(I,K) 

2089 

$  t  r  5  4  ♦  T  P 

2050 

Sr>  =  3'3  +  TMF 

2051 

5C=SC+TEMP 

2052 

SA4=3AA»TP*TP 

2053 

4 

S4B=SA9+TP*TMP 

2058 

SAC=SAC+TP«TEMP 

2055 

GO  TO  2210 

2056 

2200 

ITMp=ITMP-l 

2057 

2210 

CONTINUE 

2058 

IF  (ITFP.LT. 3)  GO  TO  2230 

2059 

C 

LINEAR  regression,  sto  dev  and  SKEW  VS  KEAN 

20b0 

TP=ITMP 

2061 

SAA=S4A-SA«SA/TP 

2062 

SAB=SA8-3A*SB/TP 

20b3 

SAC=SAC-SA*SC/TP 

2068 

C 

LIBIT  REGRESSION  COEFFICIENT  FOR  CONSISTENCY 

2065 

B8=SA0/3AA 

2066 

IF(B9.GT..25)BB=,25 

2067 

IF(P9.LT.(-.25))68»-.25 

2068 

BC=SAC/3AA 

2069 

IF  CBC.GT  ,1.)BC=1 , 

2070 

IFC8C.lt.  («1.))8C=-1. 

2071 

C 

regression  constants 

2072 

SASSA/TP 

2073 

SH=SP/TP 

2078 

CU=sa-33»3A 

2075 

SC=SC/TP 

2076 

CC=3C-BC*SA 

2077 

C 

COMPUTE  SMOOTHED  3TATI3TIC3 

2078 

00  2220  1  =  1, ND UR 

2079 

IF  (NLGG(I,K).LT.3)  GO  TO  2220 

2080 

TCMPaAvCI,K)-ALOG(P(I)) 

2081 

S0(I,<)=CB*88«TEMP 

2082 

IF  (SD(I.K).lT.O.)  S0(I,K)=0. 

2083 

SKEw(I,K)=CC*F<C*7EMP 

2088 

2220 

CONTINUE 

2085 

2230 

I p  CISKEw.LE.O)  GO  TO  2250 

2096 

00  2280  I  =  1,nOiir 

2087 

3KF.wn,K)  =  SKwCI) 

2088 

2240 

CONTINUE 

2089 

2250 

CONTINUE 

2090 

IF  (NUUR.LT.3.AN0.ISKEW.LE.0J  GO  TO  2290 

2091 

WRITE(6,50) 

2092 

WRITE(6,2260) 

2093 

2260 

for;AT(//29H  ADOPTED  FREQUENCY  STATISTICS) 

2098 

WRITE (6,890) (AA(I) ,A8(I) ,I=1,N0UR) 

2095 

00  2230  K=NSTXX,NSTA 

2096 

WRITE (6, 1070) I3TA (K) , (Av (I,K) , I  =  1,N0UR) 

2097 

WRITE (6, 1080) (33(1,10 ,1  =  1 rNCUR) 

2098 

WRITE (6,1090) (SKEW (I,K) ,1=1, NDUR) 

2099 

W,RrTE(6, 1100)  (00(1, K), 1=1, NOUR) 

2100 

IF(IPCMS.GT.O)WRTTE(7,2270) (ISTA(K),AA(I),ABCI),Av(I»K)»S0(I,K)»3K 

2  iOl 

1Ew(I,K) ,OQ(I,K) ,Anvr(I,K) ,T=1,N0UR) 

2102 

2270 

F,1HM4T(I8,1X,A3,A9,3F8.3,2F8.2/  ( 1 8 , 2 48 , 3F8 . 3 , 2F8 .2  )) 

2103 

2280 

CONTINUE 

2108 

2105 

C  * 

2106 

2290 

TMPA=100. 

2107 

X(l)=3.73 

2108 

X(2)=3.09 

2109 

X(T)=2.33 

2110 

X  £  a )  =  1 .68 

2111 

X (5) =1 .28 

22  13 
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EXHIBIT  6 


X(fe}=.52 

2iii  .  _ 

4RITE  (6,50) 

21 14 

HP  I TE  (6 » 2300) 

2115 

2300 

fqph4T(26h  computed  frequency  curves) 

2116 

On  gnro  KsNSTXX,  *'STA 

2117 

T’1PB  =  0. 

2118 

T 1-1 P  P  =  0  . 

2119 

DO  2400  11=1, NOUR 

2120 

IsNDUS-I !♦! 

2121 

IF(NL0GCI,K).LT.3)  GO  TO  2310 

2122 

THPP=TMPP+1. 

2123 

TP=SKEw(I,K) 

2124 

TflPB  =  TMP8  +  ANYR(I,K) 

2125 

2310 

DO  2390  J=1 , 13 

2126 

IF  (NLOG(I,K) .LT, 3. AND. NSTAT.LT. 1)  GO 

TO  2380 

2127 

TEMP=0. 

2120 

IF(J-7)2320, 2340, 2330 

2129 

2320 

TEMP=x (J) 

2130 

GO  TO  23«0 

2131 

2330 

TEMP=-X(14-J) 

2132 

C 

PEARSON  TYPE  III  TRANSFORM 

2133 

23ao 

IFCTP.EG.O.)  GO  TO  2370 

2134 

TEMP  =  2./'TP«C(TP/6.*(TEMP-TP/6.)  +  l.)**3 

-1.) 

2135 

TMP=(-2.)/TP 

2136 

IF(TP)  2350,2370,2360 

2137 

2350 

IF(TEMP.GT.TMP)  TEMP=TMP 

2138 

GO  TO  2370 

2139 

2360 

TT  f  TF pp , L  T . TMP )  TFNI’  =  t!,P 

2370 

TMPsAyf I,K)*TENP«SO (I,K) 

2141 

QR(J, I)=10.»*TMP-D3(I.K) 

2142 

IF(3RCJ» tJ.LT.O.)  OR  (  J » I ) 30 . 

2143 

IFdl.EQ.l.OR.J.LE.S)  GQ  TO  2390 

2144 

T"P=OR(J, I+1)»P(I)/P(I*1) 

2145 

IF(ORCJ,I).LT.TMP)OR(J,n=TMP 

2146 

GO  TO  2390 

2147 

'.300 

3RCJ, Ijs-l. 

2140 

i'39P 

CONTINUE 

2149 

2400 

CONTINUE 

2150 

IF  (TNPP.LE.O.)  GO  TO  2450 

2151 

PLTT(11=.01 

2152 

PLTT  C2)  =  , 1 

2153 

PLTTC3)=1. 

2154 

PLTT(4I=5. 

2155 

PLTT(5)=10. 

2156 

PLTTC6)=30. 

2157 

PLTT(7)=50. 

2158 

PLTTC3)=TNP4-PLTT(6) 

2159 

P  L  T  T  (9)=T4PA«PLTT(5) 

2160 

PLTT(10)  =  T’!  =  A-PLTT(4) 

2161 

PLTT(11)=Tmpa-PLTT(3) 

2162 

PLTTC2)  =  TMPA-PLTT(2) 

2163 

PLTT(1!)=TMPA-PLTT(1) 

2164 

C 

PLOT  VALUES  EXCEEOING  13  ARE 

EXPECTED  PROBABILITY 

2165 

T'lPsTNPB/TMPp 

2166 

PLTT(14)  =  .Ol*U.  +  16O0,/TMP*»1.72) 

2167 

PLT  T ( 15) =  .  1»(1. ♦200. /TKP** 1,55) 

2168 

PLTT(16j=  I.*(t.*26./TFP**1.16) 

2169 

PLTT(17)=  5.*(1.*6./TMP»*1.04) 

2170 

PLTTf lfl)310.«(l .*3./TMP**1.04) 

2171 

PLTT(19)=30.*(1.*,4(,/TMP**,925) 

2172 

PLt  T (20) =50 , 

2173 

FLTT(21)=TMPA*PLTT(19) 

2174 

plTT(22)4TmpA-PLTTC10) 

2175 

PLTT(23)=TMPA-PLTT(17) 

2176 

TLTT  (2a)=T4PA-PLTT(16) 

2177 

PLTT(25)  =  T4PA-PLTTU5) 

2170 

PLT1  (2fc)sTNPA-PLTT(14) 

2179 

4RJ7E{S,24I0)ISTA(K) 

2160 

2«10 

Frjn-.AT(/f.H  STATION  10) 

2161 

KP  I  re  (6, 2420) (AA(I),AR(I),I  =  1,NDUR) 

2162 
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2C20  FCrm.sT(4X,  16H  PLOT  EXP  PRC8  4X  ,  A3,  A4 ,  9  (  3X,  2  AO  )  ) 

O'.)  2«40  J=1 , 1 3 

WRITE (6,2430)  PLTT ( J) , PUTT ( J* 13) , (QR ( J, I ) , I»1 , NOUR ) 

2430  FORMAT (2Fl0.2,9Fll,0) 

2440  CONTINUE 
2450  CONTINUE 

N3TAT=N3TAT-NSTA 

N3TA=NSTAT 

IP(NSTAT.GT.IO)  NSTAsiO 
IP(NSTA.GT.O)  GO  TO  2140 
GO  TO  100 
END 

SUBROUTINE  CrOijT(RX) 

DIMENSION  B(10),R(10,11),RX(10,11) 
common  dtrmc, nindp, b 

NVAR3NINDP+1 
OfJ  20  J  =  l,  NINDP 
00  10  X=i,NVAR 
10  R(J,K)=RX(J,K) 

20  CONTINUE 

IF(mIN0P.GT.1)GQ  TO  30 

B(1):R(1.2)/S(l,l) 

DTRMC=B(l)*8(i) 

RETURN 

c  «**•*«  DERIVED  MATRIX  ********************* 

30  DO  40  K-2, NVAR 
40  R(1,K)=R(1,K)/R(1,1) 

DO  00  K=2, NINDP 
ITPsK-1 

DO  60  J=K, NINDP 
DO  50  1=1, ITP 
l  =  K-I 

50  R(J,N)=«tJ,K)-R(J,L)*R(L,K) 

IF(J.EQ.K)  GO  TG  60 
R(K, J)=RCJ,K)/R(K,K) 

60  CONTINUE 

DO  70  1  =  1, ITP 
LsK-I 

70  R(K, NVAR) =R(K, NY AR)-B(L, NVAR) *H(K,L) 

00  R(X,NVAH)=R(K,NVAR)/R(K,K) 

C  ******  8  aCk  SOLUTION  ********************** 

B(NINDP)=R(NINDP,NVAR) 

DO  100  1=2, NINDP 
J=NV»R-I 
I X= I -1 

n ( J ) =S ( J , NVAR) 

D -I  ~-J  L  =  1 , 1 X 
M  =  J  +  L 

90  D(J)=3(J)-8(K)*R(J,K) 

100  CONTINUE 
DT  RmC  =  0 , 

DO  110  J=l, NINDP 
110  DTS,“C  =  0TRMC*BCJ)*RX(J,NVAR) 

RETURN 

END 

FUNCTION  RNGEN(IX) 

RANDOM  NUMBER  SUBROUTINE  FOR  a  binary  machine 
GENERATES  UNIFORM  RANDOM  NUMBERS  IN  THE  INTERVAL  0  TO  1 
GENERAL  USAGE  IS  AS  FOLLOWS 
ASRNGEN(IX) 

ix  should  es  Initialized  to  zero  in  the  program 

I AHG  CAN  BE  ANY  LARGE,  ODD  INTEGER 

constants  must  ue  computed  8 y  following  eruations  .  .. 

*  *  *  *  icoNi=(2**c(an)/2))*3  *  *  *  * 

*  *  *  *  IC0N2= (2 « *B ) -1  *  *  *  * 

*  *  *  *  FC0N3=1./(2,**B)  *  *  *  * 

WHERE  B=  NUMBER  OF  BITS  IN  THE  INTEGER  WORO 

OATa  IARG/7 50021/ 

IF(IASr..EO.IX)  00  TO  10 


2164 

2164 

2105 

2106 

2107 

2108 
2109 

2190 

2191 

2192 

2193 
2)94 

2195 

2196 

2197 
2190 

2199 

2200 
2201 
2202 

2203 

2204 

2205 

2206 
2207 
2206 

2209 

2210 
2211 
2212 

2213 

2214 

2215 

2216 

2217 

2218 

2219 

2220 
2221 
2222 

2223 

2224 

2225 

2226 
222  7 
2220 

2229 

2230 

2231 

2232 

2233 

2234 

2235 

2236 

2237 

2238 

2239 

2240 

2241 

2242 

2243 

2244 

2245 

2246 

2247 
2240 

2249 

2250 

2251 

2252 
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1Y=IX 

IC0Nl=167772l9 
10  IY"IV*IC0NI 

i^0N::-2oi47  4:»7l>>'lV/o55 
IFfI-.LT.05  IYe:Y+IC0N2*l 

R -JGE N=  I  Y 

FC0N3=. 35527 t3678£«la 
K'ir,EN  =  KNGEN*?C0N3 
SE7USN 

End 


2254 

2255 

2256 
22  j  7 

2258 

2259 
22aC 
2261 
2262 
2263 
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EXHIBIT  7 
INPUT  UmTA 


A  Three  title  cards,  first  must  have  an  A  in  column  1 
B  Specification  card 

1.  NDUR  -  Number  of  durations,  dimensioned  for  8. 

2.  IYRA  -  Earliest  year  of  record  at  any  station,  dimensional 

for  100  years  (NYRS)  and  NYRS  times  NDUR  (B-l) 
dimensioned  for  4 00. 

3.  ISKEW  -  Indicator,  positive  value  calls  for  reading 

skew  coefficients  for  region. 

U.  KEEP  -  Number  of  stations  to  keep  from  the  immediately 
precedip^  job,  dimensioned  for  10. 

5.  ICONV  -  Indicator,  positive  value  calls  for  reading 

factors  to  convert  volume s  to  average  flow  rates. 

6.  IPCHQ  -  Indicator,  positive  value  calls  for  punching 

recorded  and  reconstituted  flows  on  cards. 

7.  IPCHS  -  Indicator,  positive  value  calls  for  punching 

statistics  on  cards. 

8.  NSTAT  -  Number  of  stations  for  which  statistics  axe 

to  be  read  in,  leave  blank  if  statistics  are  to 
be  computed,  no  limit  on  number. 

9.  NSMTH  -  Indicator,  blank  or  positive  value  causes 

smoothing  of  statistics. 

10.  INCAD  -  Indicator,  positive  value  calls  for  adjustment  of 
increment  to  reduce  skew  coefficient.  DO  NOT 
use  routinely  as  frequency  curves  will  be  biased. 

C  Duration  description  card 

1.  AAAB  -  Title  of  duration  such  as  "PEAK"  or  "l-DAY," 
NDUR(Bl)  items 

D  Skew  coefficients,  omit  if  ISKEW  (B3)  is  not  positive 

1.  SKW  -  Regional  skew  coefficient  for  each  successive 
duration,  NDUR(Bl)  items 

E  Stations  kept,  omit  if  KEEP(BU)  is  not  positive 

1.  KEPT  -  Station  number  (ISTA)  of  station  in  preceeding 
job,  KEEP(b4)  items.  Should  be  listed  in 
same  order  as  appearing  in  previous  job. 


EXHIBIT  7 


F  Conversion  factor,  omit  if  IC0NV(B5)  is  not  positive 

1.  P  -  Factor  by  which  flows  for  each  successive 
duration  are  divided  to  convert  to  average 
rate  of  flow,  NDUR(Bl)  items 

G  Data  cards,  omit  if  NSTAT(B8)  is  positive 

1.  ISTA  -  Station  number,  limited  to  five  digits 

2.  IYB  -  Year  number 

3.  QM  -  Flow,  NDUR(Bl)  items,  -1  indicates  missing 

record.  If  record  for  entire  year  is  missing, 
omit  card  for  that  year. 

H  Card  blank  after  Col  1  to  indicate  end  of  flow  data,  omit 
if  NSTAT(B8)  is  positive. 

I  Input  statistics,  omit  if  NSTAT(B8)  is  not  positive. 

Supply  NDUR(Bl)  cards  for  each  station  and  data  for  NSTAT 
(B8)  stations.  The  order  of  the  durations  must  be  maintained 
for  all  stations. 

1.  ISTA  -  Station  number,  limited  to  five  digits. 

2.  AAAB  -  Title  of  duration  (see  C  card.) 

3.  AV  -  Mean  logarithm  for  given  station  and  duratioi 

I4.  SD  -  Standard  deviation  of  logarithms. 

5.  SKEW  -  Skew  coefficient  of  logarithms. 

6.  DQ  -  Increment  added  to  flows  before  statistics 

were  computed. 

7.  ANYR  -  Number  of  years  of  equivalent  record. 

Four  blank  cards  with  A  in  Col  1  of  the  first  after  the 
last  job  will  cause  a  normal  stop. 
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SUMMARY  OF  REQUIRED  CARDS 
723-X6-L7350 


NSTAT 
of  NDUR 


Notes 

(1)  Omit  if  ISKEW  (B3)  is  not  positive. 

(2)  Omit  if  KEEP  (B^)  is  not  positive. 

(3)  Omit  if  ICONV  (B5)  is  not  positive. 
Omit  if  NSTAT  (b8)  is  positive. 

(5)  Omit  if  NSTAT  (B8)  is  not  positive. 


EXHIBIT  8 


